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ABSTRACT 

We investigate a theoretical magnetohydrodynamic (MHD) disc problem involv- 
ing a composite disc system of gravitationally coupled stellar and gaseous discs with a 
coplanar magnetic field in the presence of an axisymmetric dark matter halo. The two 
discs arc expediently approximated as razor-thin, with a barotropic equation of state, 
a power-law surface mass density, a ring-like magnetic field, and a power-law rota- 
tion curve in radius r. By imposing the scale-free condition, we construct analytically 
stationary global MHD perturbation configurations for both aligned and logarithmic 
spiral patterns using our composite MHD disc model. MHD perturbation configura- 
tions in a composite system of partial discs in the presence of an axisymmetric dark 
matter halo are also considered. Our study generalizes the previous analyses of Lou 
& Shen and Shen & Lou on the unmagnetized composite system of two gravitation- 
ally coupled isothermal and scale-free discs, of Lou and Shen et al. on the cases of 
a single coplanarly magnetized isothermal and scale-free disc, and of Lou & Zou on 
magnetized two coupled singular isothermal discs. We derive analytically the station- 
ary MHD dispersion relations for both aligned and unaligned perturbation structures 
and analyze the corresponding phase relationships between surface mass densities and 
the magnetic field. Compared with earlier results, we obtain three solution branches 
corresponding to super fast MHD density waves (sFMDWs), fast MHD density waves 
(FMDWs) and slow MHD density waves (SMDWs), respectively. We examine the 
m = cases for both aligned and unaligned MHD perturbations. By evaluating the 
unaligned m — case, we determine the marginal stability curves where the two 
unstable regimes corresponding to Jeans collapse instability and ring fragmentation 
instability are identified. We find that the aligned m = case is simply the limit 
of the unaligned m = case with the radial wavenumber £ — ► (i.e., the breath- 
ing mode) which does not merely represent a rescaling of the equilibrium state. We 
further show that a composite system of partial discs behaves much differently from 
a composite system of full discs in certain aspects. We provide numerical examples 
by varying dimensionless parameters (3 (rotation velocity index), r\ (ratio of effective 
sound speed of the two discs), 6 (ratio of surface mass density of the two discs), q 
(a measure of coplanar magnetic field strength), F (gravity potential ratio), £ (radial 
wavenumber). Our formalism provides a useful theoretical framework in the study of 
stationary global perturbation configurations for MHD disc galaxies with bars, spirals 
and barred spirals. 

Key words: MHD waves — ISM: magnetic fields — galaxies: kinematics and dynam- 
ics — galaxies: spiral — star: formation — galaxies: structure. 
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various complementary aspects. |Lin fc Shul dl964lll966l) and 

their co-workers pioneered the classic density wave theory 
and achieved a great success in understanding the dynami- 
cal nature of spiral galax ies dLinlll987l : iBinnev fc Tremaind 
ll987l:lBertin fc Lir]|l996T) . The basic idea of analyzing such 
a large-scale density wave problem is to treat coplanar 
perturbations in a background axisymmetric rotating disc; 
this procedure has been proven to be powerful in prob- 
ing the galactic dynamics to various extents. In such a 
model development, perturbations (either linear or non- 
linear) are introduced onto a background equilibrium to 
form local or large-scale structures and to perform stabil- 
ity analysis under various situations. Broadly speaking, ax- 
isymmetric and non-axisymmetric perturbations may lead 
to aligned configurations while non-axisymmetric perturba- 
tions can naturally produce spiral wave patterns. Theoreti- 
call y, a disc system may be treated as razor-thin for simplic- 
ity dSver fc Tremaindll996t IShen fc Loull2004bl: fLou fc Zoul 
l2004l:IShen et al.ll2005l) . with a mass density p = E(r, 6)&{z) 
under the situation where the thickness of a galactic disc is 
sufficiently small compared with its radial size scale; in this 
manner, the model problem reduces to a two-dimensional 
one. 

In theoretical model investigations, two classes of disc 
models are frequently encountered. One class is the so- 
called singular isothermal discs (SIDs), which bears a flat 
rotation curve and a constant 'temperature' wit h a diverg- 
ing su rface mass density towards the centre |e.g ., IShu et all 
( 2000)]. Since the early study bv iMestell (1963) more than 
four decades ago, this idealized theoretical SID model 
has attracted considerable interest among astrophysicist s 
in vario u s con te xts of disc d y namic s fe.g.. IZa nd dl9 76l): 
iToomrel dl977l): iLemos et alJ dl99ll) : [ Goodman fc Evans! 
dl999l): IShu et alJ d2000l); iLoul 120021) : lLou fc Fan! |2002 [) : 
Lou fc Sherd (I2003T ) : lLou fc Zoul j20o3> : lLou fc Wul 12003) : 
( 2006 1)1. The othe r class is the so called s cale- fre e 
Sver fc Tremaind dl996Tl ; IShen fc Loul d2004bl) : 
IShen et alJ d2005l) : IWu fc Loul (120061) ] and is the main focus 
of this paper. Being of a more general form for differentially 
rotating discs, a scale-free disc has a rotation curve in the 
form of v oc (the case of /3 = corresponds to a flat 
rotation curve) and a barotropic equation of state in the 
form of II = K E" where LT is the vertically integrated two- 
dimensional pressure. There is no characteristic spatial scale 
and all quantities in the disc system vary as powers of ra- 
dius r. In our analysis, the rotation velocity index f3 satisfies 
— 1/4 < f3 < 1/2 for warm discs. Furthermore, it is possible 
to construct sta tionary (i.e. lo = 0) den s ity wave patt erns in 
scale-free discs dSver fc Tremaindll996l: IShen fc LoulEoblbl 
IShen et alJ 120051: IWu fc Loull2006t) . clearly indicating that 
the pattern speed of a density wave is in the opposite sense 
of the disc rotation speed. 

Because of the self-gravity, one important technical as- 
pect in dealing with thin disc galaxies involves the Poisson 
integral, relating the mass density distribution and grav- 
itational potential. In the early development of the den- 
sity wave theory, this Poisson integral is evaluated by an 
asymptotic analysis valid in the large wavenumber regime 
of the Wentzel-Kramers-Brillouin- Jeffreys (WKBJ) approx- 
imation. Based on this, the pursuit of analytical solu- 
tions leads to a further exploration of the problem. Sev- 
eral potential-density pairs are generalized in Chapter 2 
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discs [e.e 



of IBinnev fc Tremaind (119871) . One special utility is the 
logarithmic sp iral potential-density pair first derived by 
iKalnaisI (Il97ll) . This is a powerful tool in analyzing non- 
axisymmetric pert urbati o ns in an axisymmetric background 
di sc. Furthermore. [ Qianl dl992l) made use of the techniques 
of lLvnden-Bell ( 1989) and found a larger family of potential- 
density pairs in terms of the generalized hypergeometric 
functions. Using these results for scale-free discs under our 
consideration, the Poisson integral is so lved analytically for 
both aligned and spiral perturba tions llShen fc Loull2004bt 
IShen et alJl2005HWu fc Loull2006ri . 

Magnetic field is an important ingredient in various as- 
trophysical disc systems. In the interstellar medium (ISM), 
the partially or fully ionized gases make the cosmic space an 
ideal place for applying the magnetohydrodynamics (MHD) 
equations. In many astrophysical situations involving large- 
scale dynamics, the magnetic field can be effectively con- 
sidered as completely frozen into the gas. While in some 
cases the magnetic force is relatively weak and has lit- 
tle impact on the dynamical process. There do exist cer- 
tain situations when the magnetic field plays a crucial 
role in both dynamics and diagnostics , such as in spi- 
ral ga l axies and accre t ion processe s [e.g. iB albus & Hawlev 
"19981 ) : iBalbusI d2003l): iFari fc Loul dl996l. Il997l): IShu et al 



l2000l):Eou fc Fanl |1998eJ |: ILouI d2002h : lLou fc Far] J2003D 
lLou fc Zoul (12004 l2006h : IShen et alJ 



2005) ]. In terms of MHD model developm ent, we nee d to 
prescribe the magnetic field geometry. IShu fc Lil lll997l) in- 
troduced a model in which a disc is 'isopedically' magnetized 
such that the mass-to-flux ratio remains spatially uniform 
and the e ffect of magnetic field is subsumed into two pa- 
rameters llSnu fc Lilll997l : IShu et alJfcOOd: lLou fc Wull2005l: 



IWu fc LouH2006l) . In fact. lLou fc Wul d2005l) have shown ex- 
plicitly that a constant mass-to-flux ratio is a natural conse- 
quence of the frozen-in condition from the ideal MHD equa- 
tions. In parallel, th e coplanar magnetic field also serves as 
an in teresting model dLoul2002l:lLou fc ZovtoollShen et alJ 
2005) with the magnetic field being azimuthally embedded 
into the disc system. In particular, with the inclusion of an 
azimuthal magnetic field, one can construct the so-called 
fast and slow MHD density waves (FMDWs and SMD Ws) 
dFan fc Loulll996l Il997l: lLou fc Fan|[l998al : iLolliool and 
the interlaced optical and magnetic spiral arms in the nearby 
spiral galaxy NGC 6946 are sensibly explained along this 
line. 

The stability of such MHD discs is yet another lively 
debated issue. In the WKBJ or ti ght-winding regime, the 
well-k nown Q parameter criterion dSafronovlll960] : iToomrej 
Il964l) was suggested to determine the galactic local ax- 
isymmetric stability in the absence of magnetic field. 
Meanwhile, there have been numerous studies co ncern - 
ing the global stabili t y of the disc prob l em [e.g.. pan 
jl976l): lToomrd dl977l): ILemos et al l (Il99ll) : lEvans fc Reac 
ll998dlbl)TGoodman fc Evans! <ll999l ): IShu et alJ ll2000l) l. For 
example, stationary perturbations of the zero- frequency neu- 
tral modes are emphasized as the marginal instability modes 
in scale-free discs. In this context, ISver fc Tremaina dl996l) 
made a breakthrough in obtaining semi-analytic solutions 
for station ary pertur b ation configurations in a class of scale- 
free discs. IShu et all d2000l) analyzed the stability of isope- 
dically magnetized SIDs and derived stationary perturba- 
tion solutions. They interpreted these aligned and unaligned 
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configurations as onsets of bar- like instabilities. ILoI feoosft 
performed a coplanar MHD perturbation analysis on az- 
imuthally magnetized SIDs from a perspective of stationary 
FMDWs and SMDWs. 

Our analysis on two-dimensional coplanar MHD pertur- 
bations has avoided at least two major issues. If perturbation 
velocity and magnetic field components perpendicular to the 
disc plane are allowed, it is then possible to describe Alfvenic 
fluctuations and model disc warping process. If one further 
takes into account of vertical variations across the disc (i.e., 
the disc thickness is not negligible), magneto-rotational in- 
stabilities can develop (e.g., Balbus 2003). These two as- 
pects are important and should bear physical consequences 
in modelling disc galaxies. 

In a typical disc galaxy system, the basic component 
involves star s, gases, dusts, cosmic ra ys, and a massive dark 
matter halo (ILou fc Fan! 1998a. 2003J). In terms of theoreti- 
cal analysis, it would be a great challenge to include all these 
factors into one single model consideration. While limited in 
certain aspects, it remains sufficiently challenging and inter- 
esting to consider a composite system consisting of a copla- 
nar ly magnetized gas disc, a stellar disc as well as an axisym- 
metric background of a massive dark matter halo. A seminal 
analysis concerning a c omposite disc system dates back to 
iLin fc Shu ( 19661. 11968TI who combined a stellar distribution 
function and a gas fluid description to derive a local disper- 
sion relation for galactic spiral density waves in the WKBJ 
approximation. Since then, there have been extensive the- 
oretical studies on perturbation configurati ons and stabil- 
ity pro blems of the composite disc system. Ijog fc Solomonl 
l|l984allbl) investigated the growth of local axisymmetric per- 
turbations in a compos ite stellar and gaseous disc system. 
iBertin fc Romeol 1^)88) considered the spiral modes co n- 
taining gas in a two-fluid model. IVandervoord (jlflfllal lb^l 
studied the effect of interstell ar gas on oscil lations and the 
stability of spheroidal galaxies. iRomed <ll992T) considered the 
stability of two-component fluid discs with finite thickness. 
Two-fluid approach was adopt ed into modal ana lysis mor- 
phologies of spiral galaxies bv lLowe et alJ il994f) . support- 
ing the notion that spiral struct ures are long-l asting and 
slowly evolving. lElmegreenl 1)19951) and l.Tod l|l996F) suggested 
an effective Q e ff ( or Qs-g) parameter criterion ISafronovl 
Il960l : lToomrel ll964'l for local axisymmetric two-fluid insta- 
bilities of a disc galaxy. ILou fc Fanl il998tJ> explored basic 
properties of open and tight spiral density-waves modes in 
a two-fluid model to describe a co mposite system of co upled 
stellar and gaseous discs. Recentlv. lLou fc Shenl 120031) stud- 
ied a composite SID system to derive stationary global per- 
turbation configurations and further explored the axisym- 
metric instability properties IShen fc Loull2003h where they 
proposed a fairly straightforward D criterion for the ax- 
isymmetric instabilit y problem for a composite SID system. 
IShen fc Loul l2004bl) extended these analysis to a compos- 
ite system of two scale-free discs and carried out analytical 
analysis on both aligned and logarithmic spiral perturba- 
tion configurations. By adding a coplanar magnet ic field to 
the background composite SIDs. ILou fc Zoul 120041) obtained 
MHD perturbation configurations and further stud ied the 
axisymmetric instability problem llLou fc Zoj2006Tl . 

The main objective of this paper is to construct global 
scale-free stationary configurations in a two-fluid gaseous 
and stellar disc system with an embedded coplanar mag- 



netic field in the gas disc. Meanwhile, an axisymmetric in- 
stability analysis is also performed. Ther e are several new 
features compared with previous works l)Lou fc Zoull2004l : 
IShen fc Lor3l2004bt IShen et al.ll2005l) that may provide cer- 
tain new clues to understand large-scale structures of disc 
galaxies. This paper is structured as follows. In Section 2, we 
present the theoretical formalism of the problem; both the 
stationary equilibrium state and the linearized MHD pertur- 
bation equations are summarized. In Section 3, we perform 
numerical calculations for the aligned perturbation configu- 
rations. Both the dispersion relation and phase relationship 
between density and the magnetic field perturbations are 
deduced and evaluated. In Section 4, we apply the same 
procedure to the analysis of global logarithmic spiral con- 
figurations. The m = marginal stability is also discussed. 
Finally, we summarize and discuss our results in Section 5. 
Several technical details are included in Appendices A, B, 
and C. 



2 FLUID-MAGNETOFLUID DISCS 

We adopt the fluid-magnetofluid formalism in this paper to 
construct large-scale stationary aligned and unaligned copla- 
nar MHD disturbances in a background MHD rotational 
equilibrium of axisymmetry |Lgj^fc^hen| 12223!) ;|Lgj^&^23 
|2004 l: IShen fc Loul |2004blk ILou fc Wul 120051] i: ILou fc Zoul 
feoOffl : IWufc LoiJ |2006h ]. All the background physical 
quantities are assumed to be axisymmetric and to scale as 
power laws in radius r. Specifically, the rotation curves bear 
an index of —j3 (viz., v oc r _/3 ) and the vertically integrated 
mass density has the form Eo oc r~ a with —a being another 
index. Physically, the magnetofluid formalism is directly ap- 
plicable to the magnetized gas disc, while the fluid formal- 
ism is only an approximation when applied to the stellar 
disc, where a distribution function approach would give a 
more comprehensive description. For our purpose of mod- 
elling large-scale stationary MHD perturbation structures 
and for mathematical simplicity, as well as the similarit y 
between the two sorts of descriptions IShen fc Loull2004bl) . 
it suffices to work with the fluid-magnetofluid formalism. 

In this section, we present basic MHD equations of the 
fluid-magnetofluid description for a composite system con- 
sisting of a scale-free stellar disc and a coplanarly magne- 
tized gas disc. In our approach, the two gravitationally cou- 
pled discs are treated using the razor-thin approximation 
(i.e., we use vertically integrated fluid-magnetofluid equa- 
tions and neglect vertical derivatives of physical variables 
along z) and the two-dimensional barotropic equation of 
state to construct global MHD perturbation structures. A 
coplanar magnetic field is involved in the dynamics of the 
thin gas disc following the basic MHD equations. The back- 
ground state of axisymmetric rotational equilibrium is first 
derived. We then superpose coplanar MHD perturbations 
onto the equilibrium state and obtain linearized equations 
for MHD perturbations in the composite disc system. 

2.1 Ideal Nonlinear MHD Equations 

By our conventions, we use either subscript or superscript 's' 
and 'g' to denote physical variables in association with the 
stellar disc and the magnetized gas disc, respectively. For 
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large-scale MHD perturbations of our interest at this stage, 
all diffusive effects such as viscosity, resistivity, thermal con- 
duction, and radiative losses etc. are ignored for simplicity. 
In cylindrical coordinates (r, 6, z) and coincident with the 
2=0 plane, we readily write down the basic nonlinear ideal 
MHD equations for the composite disc system, namely 
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for the magnetized gas disc, and 
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for the stellar disc in the 'fluid' approximation, respectively. 
Here, E denotes the surface mass density, u is the radial 
component of the bulk flow velocity, j = rv is the spe- 
cific angular momentum in the vertical z-direction and v 
is the azimuthal bulk flow velocity, n is the vertically inte- 
grated two-dimensional pressure, B r and Bg are the radial 
and azimuthal components of the coplanar magnetic field 
B = (B r , Be, 0), and (j> is the total gravitational potential. 
This cj> can be expressed in terms of the Poisson integral 



F<f>(r, 9) = j dkp 



-GE(r', ip, t)r'dr' 



[ r /2 + r 2 _ 2r 'r cos(v? - 8)} 1 / 2 



(7) 



where G is the gravitational constant, E = E 9 + E s de- 
notes the total surface mass density, and F is defined as 
the ratio of the gravitational potential arising from the 
composite discs to that arising from the entire system in- 
cluding an axisymmetric massive dark matter halo that is 
presumed not to respo nd to the coplana r MHD perturba- 
tions i n the disc plane |Syer fc Tremain5 <|l996f); |Shu et"aL 
|2000D.lLou fc Shenl j2003T l : lLou fc Zoul J20041) : IShen fc Lou 
l|2004bD : Lou & Wu (2006); Lou & Zou (2006)] with F = 1 
for a full composite disc system and < F < 1 for a com- 
posite system of partial discs. 

For a coplanar magnetic field B = (B r , Bg, 0) in cylin- 
drical coordinates (r,9,z), the divergence-free condition is 



d{rB r ) dBg 
dr d9 



= 



The radial and azimuthal components of the magnetic in- 
duction equation are 



dB r Id 

-dT = rde {u Be 



V 3 B r ) 



V 9 B r ) 



(9) 



(10) 



Among JHJ — 1101 . only two of them are independent. The 
barotropic equation of state for the scale-free discs is 



n 1 = KiY? 



(11) 



where Ki > is a constant proportional coefficient and 
rii > is the barotropic index with subscript i denoting 
either g or s for the two discs (we use this convention for 
simplicity). An isothermal equation of state has m = 1. For 
the stellar disc and the magnetized gas disc, Ki are allowed 
to be different, but we need to require n g = n s = n in order 
to meet the scale- free requirement (see the next subsection). 
It follows that the sound speed a.i (in the stellar disc the ve- 
locity dispersion mimics the sound speed) for either disc is 
readily defined by 



(12) 



where the subscript o denotes the background equilibrium. 



2.2 Axisymmetric Background MHD Equilibria 

For the axisymmetric rotational MHD equilibrium, we set 
u 3 = u s = and all terms involving d/dt and d/d6 in equa- 
tions — Hot to vanish. We now determine the background 
axisymmetric equilibrium with physical variables denoted 
by subscript V. In our notations, a background equilibrium 
is characterized by the following power-law radial scalings: 
both surface mass densities Eg, E 9 oc r~ a carrying the same 
exponent index a yet with different proportional coefficients. 
The disc rotation curves v$ and take the power-law form 
of cx r - ' 3 and therefore the z— component specific angular 
momentum jo = rvo oc r 1- ' 3 for both discs. In our composite 
MHD disc system, a background magnetic field is purely az- 
imuthal in the gas disc, that is, B® = and Bg = Bo oc r~"' . 
A substitution of these power-law radial scalings into equa- 
tions Q — 111 2t leads to the following radial force balances 



i 2 g + rd(j) /dr + (1 - ~f)C\ 



(13) 



in the magnetized gas disc, where C'a is the Alfven wave 
speed in a thin magnetized gas disc defined by 



1 



4^E 9 



Bgdz 



(14) 



varying with r in general, and 

vq 2 = —ota 2 . + rd(j)o/dr (15) 
in the stellar disc. Poisson integral Q yields 
F<j, = -27rGrE Po(a) , (16) 
where the numerical factor Vo(a) is explicitly defined by 
r(-a/2 + l)r(a/2 - 1/2) 



Po(a) 



2r(-a/2 + 3/2)r(a/2) 



(17) 



Coupled Magnetized Scale-Free Discs 5 



with r(- ■ ■ ) being the standard gamma function. Expression 
11711 can also be included in a more general form of Vm(/3) 
defined later in equation (|55^ , 

In order to satisfy the scale- free condition, radial force 
balances (1131 and 11511 should hold for all radii, leading to 
the following simple relation among the four indices 



2/3 = a(n - 1) = 2j - a 



(18) 



This in turn immediately gives the explicit expressions of 
indices a, 7 and n in terms of 0, namely 



a = 1 + 2(3 



7 : 



(1 + 4/3) 



(1 + 4/3) 
(1 + 2/3) 



(19) 



Since n > for warm discs, we have /3 > —1/4. Further- 
more, Poisson integral Q converges when 1 < a < 2 which 
corresponds to < (3 < 1/2. In addition, for a finite total 
gravitational force a larger /3-range of — 1/2 < (3 < 1/2 is 
allowed. Therefore, th e physical range for i s constrained 
by -1/4 < < 1/2 JSyer fc Tremainel[l996l : IShen fc Loul 
l2004bl:lLou fc Zoull2004l: IShen et alJl2005tl . 

For simplicity, we introduce the following parameters 
Ai, Di, and q defined by 



A 



cl 



(1 + 2/3)r 2 f 



2 2 
q a g = 



"0 



AiDi 



(20) 



l + 2/3)r 2 " 



where the subscript or superscript i denotes either g or s 
for either the magnetized gas disc or the stellar disc, re- 
spectively. We see that Ai is related to a scaled sound 
speed (or velocity dispersion of the stellar disc), the con- 
stant Dt — «o/[ai(l + 2/3) 1 ^ 2 ] is a scaled rotational Mach 
number, and the constant parameter q represents the ratio 
of the Al fven speed Ca to the sound speed in the magnetized 
gas disc dShen et al.ll2005l) . 

It is then straightforward to express other equilibrium 
physical variables in terms of Ai, Di and q. For the disc 
angular speed Qi and the epicyclic frequency Ki, we have 
respectively 



magnetized gas disc is typically less than the stellar velocity 
dispersion (regarded as an effective sound speed). 

With these notations, we now have from condition 1221 



Y 9 — 



So 



F^ 2 [D 2 + l + (4/3- 



l)q 2 /(4/3 + 2)]5 



2ttG(2/3Po)(1 + <5)r< 1 + 2 ' 3 ) 
fM 2 [.D 2 + l + (4/3-l)g 2 /(4/3 + 2)] 
2ttG(2/3Po)(1 + <5)r( 1 + 2 < 3 ) 



(23) 



Note t hat expressio ns 122H and 12311 reduce to expression 
(14) of lShen fe Loul (2004bj) when the magneti c field van- 
ishes ( i.e., q = 0) and are also in accordance with lShen et alJ 
(2005) when a single magnetized scale-free gas disc is con- 
sidered. 

Since the magnetic field effect is represented by the q 
parameter multiplied by a factor (1 — 7) or (4/3 — 1) [see 
cq. 113H or eq. 122H 1. we know that in the special case of 
— 1/4, the Lorentz force vanishes due to the cancella- 
tion between the magnetic p ressure and tensi on forces in 
the background equilibrium dShen et al.ll2005l) . Moreover, 
when —1/4 < < 1/4, the net Lorentz force arising from 
the azimuthal magnetic field points radially inward, while 
for 1/4 < < l/2 ; the net Lorentz force points radially 
outward iShen et alJl2005l) . 

Another point should also be noted here. From equa- 
tion 1221 . there exists another physical requirement for the 
rotational Mach number D g , i.e., it should be large enough 
to warrant a positive D 2 . This requirement for D 2 is simply 



Di 



4/3 + 2 



1 2 

q 



- 1 > 



(24) 



In the presence of magnetic field, the relation between 
D g and D s becomes somewhat involved, as c o mpared to the 
hydrodynamic case studied bv lShen fc Loul £2004b). Since 
77 1, we see from equation 1241 that when ^ 1/4, D 2 
remains always smaller than Dg, while for > 1/4, Dg is 
smaller than D 2 when 77 is of order of 1. 



A,Dj 



2£li d(r 2 Oi) 



dr 



1/2 
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A combination of equations 113H . I15L 1161 and 1201 1 leads 
to the relation among these parameters in the form of 



Di + i + 



4/3 



4/3- 



A 2 s (Dt + l) 



(22) 



= 2nG20V o T, o r 1+2l3 /F 



We now introduce two dimensionless parameters to com- 
pare the properties of the two scale-free discs. The first one 
is simply the surface mass density ratio 8 = Eq/Eq, and the 
second one is the square of the ratio of effective sound speeds 
77 = a 2 /a 2 = A 2 g /A 2 . For disc galaxies, the ratio 8 can be 
either greater (i.e., younger disc galaxies) or less (i.e., older 
disc galaxies) than 1, depending on the type and evolution 
stage of a disc galaxy. Meanwhile, the ratio rj can be gen- 
erally taken as less than 1 because the sound speed in the 



2.3 Coplanar MHD Perturbations 

In this subsection, we use subscript 1 along physical vari- 
ables to denote small coplanar MHD perturbations. For ex- 
ample, E 9 = Eq + Ef for gas surface mass density with Ef a 
small perturbation quantity under consideration. Similar no- 
tations hold for other physical variables. Now from equations 
Q — HOI and the specification of the background rotational 
MHD equilibrium, we readily obtain linearized equations for 
coplanar MHD perturbations below, namely, 



ggf 
dt 



+ 



ld_ 



(rE 9 <) + a 



9E? Egdj? 



+ 



= 



(25) 



dt 



+ Qg 



09 



r dr 



2 ^1 

^0 



(1 - 4/3)ClE? 



^0 



2Egr 
dzBo 
4nr 



d{rbg) 



dzbg d(rBo) 
47rr dr 

8br 

00 



(26) 
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y3 
! ^1 



rfzfc r d(rB ) 
4tt dr 



(27) 



in the magnetized gas disc and 



2 nA = - d 



dt 89 r dr 

at + 2«. 1+ • a# ~ ae 



vis 

2 1 . j 



(28) 



in the stellar disc, respectively. All dependent variables are 
taken to the first-order smallness with nonlinear terms ig- 
nored. The perturbed Poisson integral appears as 



dip 



-G(Sf + Y,\)r'dr' 



[ r /2 + r 2 _ 2r'r cos(<p - 0)]V2 



(29) 



Note that for a gravitational potential perturbation, we no 
longer have the factor F by the simplifying assumption that 
the dark matter halo exists as a background and does not 
respond to perturbations in the two discs. Together for a 
coplanar magnetic field perturbation b = (b r ,bg,0), we have 



d(rb r ) dbg_ _ 

dr + 80 ~ ' 
db r 13, 



(30) 



from the divergence-free condition and the magnetic induc- 
tion equation. 



2.3.1 Non-Axisymmetric Cases ofm^l 

As the background equilibrium state is stationary and ax- 
isymmetric, these perturbed physical variables can be de- 
composed in terms of Fourier harmonics with the periodic 
dependence exp(iwi — imff) where u> is the angular frequency 
and m is an integer to characterize azimuthal variations. 
More specifically, we write 

Ei = S* (r)exp(iurf — imO) , u\ = U 1 (r)exp(io;t — imff) , 
j{ — J 1 (r)exp(iuit — imff) , tf>i = V(r)exp(iut — imO) , 
b r — i?(r)exp(kjt — imO) , be = Z(r)exp(iujt — imff) , 

(31) 

where we use the italic superscript i to indicate associations 
with the two discs and the roman i for the imaginary unit. 
We also define V 1 (i = s, g) for the corresponding gravi- 
tational potentials associated with the stellar and gaseous 
discs such that V(r) = V s (r) + V 3 (r). By substituting ex- 
pressions 1311 into equations (12511 — 1301 . we readily attain 



i(«-mfi fl )5 9 + i4-(r£gI/ s )-imEg^ =0 
r dr " u r 2 



(32) 



i(w - mQ, g )U 9 - 2Q g 



J 9 _ d* 9 (1 - 4/3)Cjg g 
r ~~ dr + 2E?r 



yg 



dzBo 
Anr 



d{rZ) 
dr 



+ imR 



1 

V 9 



i(w - m£l g )J 9 + ^-U 9 = im* s + ^ 

for the magnetized gas disc, and 

i(w-rrfo a )S' + -4-{rT?oU a 
r dr 



dzZ d(rB ) 
Airr dr 

dzR d(rB ) 
4tt dr 



(33) 



(34) 



imEo- 



i(w - mSl s )U s - 2D. 



,1 s 



s — 



dm 3 

dr 



(35) 



i(w - mD s ) J s + ^-U s = im* s 

for the stellar disc, respectively, where we define 
¥ = of^VEj + V . 

The perturbed Poisson integral now becomes 
-G{S 3 +S s )cos(mip) r'dr' 



V(r) 



d(p 



10 ( r ' 2 + r ' 2 ~ 2r'r cos ip) 1 / 2 
and perturbation magnetic field variations reduce to 
d(rR) 



(36) 



(37) 



dr 



imZ — , 



(38) 



i(w - mD g )R + '-^-U 3 = , 
r 

iuiZ = 4-(rn g R- B U 3 ) . 
dr 

A combination of equations (l.'ilifi — l|35[l and I37[l — l|38fl consti- 
tutes a complete description of coplanar MHD perturbations 
in a composite disc system. From equation 1381 . we obtain 

mB U 3 



Z = 



i d(rR) 
m dr 



R 



r{to — mO, g ) 



(39) 



A substitution of expressions 1391 into the radial and az- 
imuthal components of the momentum equation (IMol and 
O yields 

J 9 _ d* 9 (1 - A(i)C\S 9 



i(w - mD g )U 9 - 2Q g ^- = + 
1 



-cl 



d*_ 

dr 2 



r dr 
12/3 d 2/3(1 + 4/3) - m 2 



2E 9 r 



2r dr 



iU 3 



LO — mfi 



i(ui — mD. g )J 3 + 



" -J- ^-U 9 



2Q a 



im* 9 



mC 2 A (1 ~4P)U 3 
2r(uj — mSl g ) 



(40) 



(41) 



We are primarily interested in stationary global configura- 
tions of zero-frequency MHD perturbations yet without the 
usual WKBJ approximation. By setting uj = 0, we can pre- 
serve the scale-free condition. Stationary perturbation pat- 
terns inourh^meofreference have been studied previously 



fe.g..lSver fc Tremaind ll996^:IShu et alJl|200(^ : lLoulj2002}):_ 

3ul I 



IShen fc Loul J2003t . l2004bl) : Ilqu fc Zouf 120041) : IShen et alJ 
i2005Tl : Lou & Wu (2005)] and our analysis here is more 
general, involving a composite system of two scale-free discs 
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with a coplanar magnetic field. For m 1, we set u> 
MHD perturbation equations to deduce 



mtlgS 3 + -^-(rE 9 iU 3 ) + mE 9 ^ 
r dr r 2 



, 



in 



(42) 



(51) 



mQgiU 3 + 2tl g 



m 



dr 2 



+ 



J 9 _ d* 9 _ (1 ~4/3)C'lS 9 
~7 ~ ~dr 2E 9 r 

1-12/3 d 2/3(1 + 4/3) - rn 2 



2r dr 



it/ 9 



mQ,„J 3 + 



20, 



■it/ 9 



-m* 9 + 



Cl(l-4/3) it/ 9 



2r 



a, 



(43) 



(44) 



"9 "9 

for the magnetized gas disc. Meanwhile, we derive from two 
independent equations of I39H 



d fBoiU 



dr \ mfl a 



B U 9 



R= ^~' (45) 

'9 / '"9 

for the magnetic field perturbation. For coplanar perturba- 
tions in the stellar disc, we derive 



mSl a S 3 + ~^(rEgi£/ s ) + mEg^ 







(46) 



r dr 

mfl s iU 3 + 2Q S — 
r 

mn s r + ^-iu 3 

As part of the derivation, we deduce from the last two equa- 
tions of H46> for the expressions of U 3 and J 3 as 

-imfi s 

r 



v 



d*r 

dr ' 

-m^ 3 



m 2 Q, 2 — k 2 V r dr I 



J s 



rn 2 Q, s k? d \ T , 
1 + ttt^-t- * 



(47) 



m 2 f2 2 _ K 2 y r 1 2J1 S dr y 



A substitution of equation l|47|l into equation l|46|l gives a 
single equation relating S 3 and ^ s for coplanar perturba- 
tions in the stellar disc, namely 



mQ. s S 3 + 



ld_ 

r dr 



mf2 s rEg (2 



mEf, 



+ 



2 d 



(48) 



* s = 



/tj \ r' 2rf2 s dr, 

Now equations (H2l — AH . lESI together with equation (f^Tt 
form the coplanar MHD perturbation equations for con- 
structing global stationary non-axisymmetric m ^ 1 con- 
figurations for both aligned and logarithmic spiral cases in 
a composite MHD scale-free disc system. 

2.3.2 The Axisymmetric Case of m — 

We now examine the axisymmetric case of m = for which 
equations 1421 — 1441 . 1481 1 become degenerate. Instead, we 
adopt a different limiting procedure by setting m — and a 
small u> =fc in equations ETl . and l!35t to 

obtain 



^S 9 = i-f(rE 9 i[/ 9 ) , 
r dr 



Lj 2 iU 3 — 2u}Q„— = — 



d* 9 (1 - 4/3) CIS 9 
"-fr*" 2E 9 r 



dr 2 



1 - 12/3 d 2/3(1 + 4/3) 

^ 7 "T" O 



2r dr 



(49) 



(50) 



if/ 9 



dr 



Boit/ 9 



(52) 



for the magnetized gas disc. For coplanar perturbations in 
the stellar disc, we have in parallel 



wS 3 = I-^-(rEgi[/ s 
r dr 



uiiU 3 

J 3 



2fi 



r 



dfr 3 
dr 



(53) 



20 s 



-iU 3 



where the gravitational coupling between the two scale-free 
discs is implicit in ty 9 and W expressions. 

For global stationary axisymmetric MHD perturba- 
tions, we t ake the limiting procedure u> — > in equations 
t — 1531 llLou fc Zoul2004l:IShen fe Loii2004ri IShen et alJ 



2005: [Lou fc Zoull2006ft rather than setting u = in MHD 
perturbation equations and then let m — » 0. Note that these 
two procedures of taking limit lead to dif f erent results except 
for certain special cases [e.g.. IShu et alJ J2000f) l . 



3 GLOBAL STATIONARY ALIGNED MHD 
PERTURBATION CONFIGURATIONS 

In the following two sections, we solve the Poisson integral 
by introducing two kinds of scale-free global MHD pertur- 
bation patterns, namely, the aligned and the logarithmic 
spiral configurations. These two kinds of MHD perturbation 
structures are the results of the special forms in a l ist of 
the p otential-density pairs feinnev fc Tremainelll987t iQianl 
1992). For the aligned case, the maximum density pertur- 
bations at different radii line up in the azimuth, while for 
the spiral case, the maximum density perturbations involve 
a systematic phase shift at different radii such that a loga- 
rithmic spiral pattern emerges. 

In this section, we derive aligned stationary perturba- 
tion dispersion relations on the basis of the results of the pre- 
vious section and discuss solution properties in detail includ- 
ing the dependence of disc rotation speed on several param- 
eters, the phase relationship among the surface mass density 
perturbation in the two discs and the perturbed magnetic 
field. As an example of illustration, we mainly focus on a 
composite system of full scale-free discs (i.e., F = 1). 



3.1 Dispersion Relation for Global 

Aligned Coplanar MHD Perturbations 

For aligned MHD perturbations, we select those perturba- 
tions that carry the same power-law variations in r as those 
of the background equilibrium, that is, in expressions 13H 



S\r) = eV" 1 - 2/3 



V i (r) = -2nGrS i V m (0) , (54) 



where e l is a sufficiently small amplitude coefficient, super- 
script i denotes either g or a for either gas or stellar discs, 
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respectively. The numerical factor V m (0) is defined explic- 
itly by 



Vm(0) 



T(m/2 -0 + l/2)r(m/2 + (3) 



(55) 



2r(m/2 -/3+ 1 )r(m/2 +0 + 1/2) ' 

wher e -m/2 < < (m + 11/2 iQianll992l:ISver fc Tremainei 
Il996l: IShen fc LoulEooibl IShen et alJl2005h . Note that the 
requirement of— 1/4 < < 1/2 automatically satisfies this 
inequality for m ^ 1. In the limit — » 0, we simply have 
P m = 1/m. For m = 0, equation 1551 reduces exactly to 
equation 11711 by setting a = 1 + 20; meanwhile, we also 
have the limiting result of 20Vq(0) — > 1 as — > 0. 

We start with equations - 6H 1 Hit an d by 
applying power-law variations that satisfy scale-free condi- 
tions. Here, we have U oc r - ' 3 and J oc r 1_/3 . By proper 
combination and simplification, we then arrive at 



m o S 9 -^^iU 9 + ^J 9 







2mO, g r 2 



2Q n 



4irG0T„ 



mflnJ 9 + 



40a 2 g + (1-40)C 2 A 



2rEg 



(56) 



(57) 



S 9 + 47vG0V m S s , 



Z=- 



(1 - 0)rQ g - 

2TYGVm - 



2 J mSl g r ' 



C%(l-40) 



2rQ a 



iU 9 = 



rEg 



5 s + 2TvGmrV m S B 



R = 



B U 9 



for the magnetized gas disc and 
r 2 (m 2 -2 + 20)tl 2 



(m 2 +10- 10 2 ) 



at + 27rGrEoP„ 
+ 2TrGrT, s VmS 9 



(58) 



(59) 



(60) 



for the stellar disc. Note th at equations 1561 — 1591 are 
the same as equation (43) of IShen et alJ (j5x)5) for a sin- 
gle magnetized scale-free disc by setting S s — and equa- 
tion 160 1 is exactly t he sam e as the first one of equation 
(35) in IShen fc Loul <2004bl) . A combination of equations 
1561 — 1581 and 16U1 then gives a complete solution for the 
stationary dispersion relation in the gravitational coupled 
MHD discs that are scale free. Because these equations are 
linear and homogeneous, to obtain non-trivial solutions for 
(S 9 , S" ,U 9 , J 9 ), the determinant of the coefficients of equa- 
tions 1561 — 1581 and 1601 should vanish. This actually gives 
rise to the stationary MHD dispersion relation. Meanwhile, 
in order to get a physical sense of the dispersion relation, we 
solve the above MHD equations directly. 

A combination of equations 1561 and 1581 produces re- 
lations of ill 9 and J 9 in terms of S 9 and S s , namely 



iU 9 = ^f[(QA + E A )S 9 + QS S ] , 

J 9 = %£{[-! + 3/3(6 A + E A )]S 9 + 30S A S S } 



(61) 



where for notational simplicity, we define 
6a = (27rGrP m Eg)/AA , 

E A = (-a 2 3 + Qlr 2 )/A A , (62) 
Aa = (1 + 20)tl 2 g r 2 + (20 - 1/2)C\ . 

Here, we use the subscript A to indicate the 'aligned case'. 
In particular, Q A and Ha are two dimensionless constant 
parameters. A substitution of equation 1611 into equation 
1571 and a further combination with equation 1601 lead to 

[/Ca(6a + Ha) + C a - 2-KGY 1 9 ) r{20T m )]S 9 

+ [K. A Q A - 27vGX 9 r(20T m )]S s = , (63) 
2TvGT,t ) rV m S 9 + [Ma + 2nGE a rP m ]S s = , 
where for notational simplicity we define 
Ka = ijn 2 + 60)Q 2 g r 2 - [m 2 - (1 - 10) 2 /2\C\ , 
La = (1/2 - 20)C 2 a - 2Vl 2 g r 2 + 20a 2 g , (64) 
Ma = (m 2 + 20- 2)Q 2 a r 2 /(m 2 +40- 40 2 ) - a 2 . 

We now obtain the stationary MHD dispersion relation by 
calculating the coefficient determinant of equation 1631 . Af- 
ter a proper rearrangement, we obtain 

^-[/Ca(6a + Ha) +Ca - 27rGEgr(2/3P m )]x 

(Ma + 2^GE^P m ) = 4^ 2 G 2 EgE^r 2 P^ (l - 

(65) 

where the left-hand side consists of two factors. One can 
show that the left bracket is exactly the dispersion relation 
for a single coplanar magnetized scale-free disc discussed 
by IShen et al.l i2005T) and the second parentheses denotes 
the dispersion relation for a single hydrodynamic scale-free 
disc. The right-hand side denotes the effect of gravitational 
coupling between perturbations in the two scale-free discs. 
By setting = for an isothermal composite disc system, 
equation (|651 the n reduces to dispersion relation (62) of 
iLou fc Zoul a2004f) where MHD perturbations in isothermal 
fluid-magnetofluid discs are investigated. 

These calculations are straightforward but tedious and 
we turn to the following subsections for further analyses. 



3.2 The Aligned m = Case 

As already discussed earlier, the m = case is special and 
should be treated in the procedure of setting m = and 
then taking the limit of uj — » 0. By letting u) — > in equa- 
tions gsj - E3> 

we find U = and R = similar to 
the earlier work [e.g . Ishu et al.l |2000D ; lLou fc ShenH2003t) : 
IShen fc Loul <l2004bfl : lLou fc ZoulfcOothlShen et alTj2005l) 1 . 

By further requiring the scale-free condition, other physical 
quantities should be in the forms of Z oc r~ 7 and J oc r 1- ' 3 
which are exactly the same as those of the background equi- 
librium. Since these perturbations are axisymmetric, it turns 
out that such perturbations simply represent a sort of rescal- 
ing of the axisymmetric background equilibrium. Neverthe- 
less, by carefully taking limits in calculations, we can also 
find a stationary 'dispersion relation'. We leave this anal- 
ysis to Appendix A for a further discussion and focus our 
attention onm) 1 cases in the next subsection. 
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3.3 The Aligned m > 1 Cases 

3.3.1 Dg Solutions for the Dispersion Relation 

We aim at evaluating the dispersion relation numerically to 
explore the dependence of Dg or Dg on a set of dimension- 
less parameters (i.e., m, /3, n, 5, q, F). One way is to begin 
directly from equation 16511 . and the other way is to cal- 
culate the coefficient determinant of equations 1561 — 15811 
and 1601 . For both, we use equations 1201 — 12311 to make 
systematic calculations. Making a full use of mathemati- 
cal tools for computations, we choose to follow the latter 
procedure. This is a straightforward but onerous task; we 
show the results below and leave the technical details to 
Appendix B. Here we first compute Dg because its expres- 
sion appears simpler, and then obtain the corresponding D 2 
using equation 1221 . When D 2 S is smaller than Dg, we show 
D 2 solutions in figures. Otherwise, we show both Dg and 
Dg in figures to identify physical solutions. We now intro- 
duce a few handy notations below to simplify mathematical 
manipulations, namely 

Ai{[3) = m 2 +4/3-4/3 2 , 

B m (fi) = (1 + 2/3) (m 2 -2 + 2/3) , 

C(J3) = 2/3P /[F(l + 2/3)] , 

Hi(P) = AiVm/C + B m . 

In the following, the stationary MHD dispersion relation is 
written as a cubic equation in terms of y = Dg, namely 

C 3 V + C?y 2 + C A y + C$ = . (67) 

The four coefficients C A ,C A , C A , and C A are functions of 
dimensionless parameters m, /3, 5, r\, and q 2 and are defined 
explicitly by 

C 3 = BmTL m rj , 



(66) 



(Br, 



A A \T-/ A I (Am + B m )(Hi 
*^m ) i^m T 



B m )5 



(1 + 5) 

(A A + B m )(H A 5 + B m ) , n A 2 

r ^2Ql 



— A A H A 4- 

•s^m ' *-m I 



(1 + 5) 
(At + B m )(Ht 



B m )5 



(1 + 5) 
(At + B m )(HiS + B m ) 



(68) 



r1 A 2 
+ C 1Q q 



(1 + 5) 

/-iA /-iA 6 | /-iA 4 | ^A 2 

C) = °0Q39 + ( - , 0Q2 ( Z + ^0Ql9 > 

where the five new relevant coefficients are further defined 
explicitly by 

m 2 (4/3 ■ 



c 



2Q 



1) - 2At + 8/3 2 -20/3 + 6 A ,*3 r\ 
{rim. + Om0)n 



2(1 + 5) 

4(1 - 2(3)A A - (1 + 2/3) BmV, 
(4/3 + 2)(l + 5) C 



-Sn 



.4^ + 5/3-2 
(1 + 5) 



(H m + BmS)n 



_ Aj (4/3 -3) +2/3+ I 
(4/3 + 2)(l + 5) 

_ .44 (4/3 -3) +2/3+ I P., 
(4/3 + 2)(l + 5) C 

+ (Ai + Bm)(Ai + 5[3-2) 



■xjA . jo r 2BmP m S S ' 
H m + OmO - | 77 



0Q3 



c, 



0Q2 



+ 



(4(3~l) 3 (H A +BmS)v 
2(4/3 + 2) 2 (f + 5) 

(4/3-l) 2 (2^4-2/3-l)g m P r , 

C 

-2/3 



(4/3 + 2) 3 (l+5) 
(4/3-1 



~ (4/3 + 2) 

x + B, 
(4/3 - 1)(2^4 



(4/3 - 1) + (2A A 



■8r] , 



(4/3 + 2)(l + 5) C 



(4f3-l) 2 (Hm.+Bm5) V 



2/3- 



(4/3 + 2) 2 (l + 5) 



2(4/3 + 2)(l + 5) 

m ( -' 



2BmVm8 
C 



0Q1 



. (2At -2/3-1) 
■ (4/3 + 2)(l + 5) 

(2At - 2/3 - 1) 



ojA . „ c BmPmS 

H m + o m 5 — I n 



C 

1 + 5-^5) 



(69) 



Br, 



(4/3 + 2)(l + 5) 

Before solving cubic equation l|67|l . we can examine qualita- 
tively several necessary requirements between y = D 2 and 
other parameters involved. For q — 0, we have C A = 0; 
meanwhile, C A , C A , and C A become exactly the same as 
C*2, Ci , and Co, respectively in equation (39) of lShen fc Loul 
(2004b). With the introduction of a coplanar magnetic field, 
a new solution (SMDWs) emerges and in the limit of q —> 0, 
this solution shrinks to zero. By letting 77 = and 5 — » 00, we 



of a single MHD scale-free disc studied bv lShen et alJ $2(1 
[see their coefficient expressions (49)]. It is the presence of 
a stellar disc in the model that results in an extra possible 
mode for stationary perturbations (Lou & Fan 1998b). 

It is also useful to express this relation in terms of D 2 
because a physical solution should have D 2 . > and D 2 > 
simultaneously. We deduce from equation 1221 that 



(70) 



D 2 = D 2 / V + X , 

*=(l-r;)/77-(4/3-l)g 2 /(4/3 + 2) . 

Substituting equation 1701 into equation (1691 . we immedi- 
ately obtain a new relation in terms of y = D 2 , namely 



/~,A' IZ 1 ,-yA' 12 , ^A' 1 . /-tA 1 n 

Qi y +C 2 y + Cx y + C = , 

where the four primed coefficients are defined by 

Of' = Ci/rf , 

Cf = (C 2 A + ZCiX)/n 2 , 

C A ' = (C? + 2C A X + 3C A X 2 )/n , 

Cq 4 ' = C5 4 + C A X + C^X 2 + C^X 3 . 



(71) 



(72) 



It would be straightforward but tedious to express the 
coefficients of cubic equation 1711 explicitly as we did be- 
fore. We would use equation 1671 for D 2 g computations and 
equation 1711 for theoretical inferences later. 



3.3.2 Dependence of D s and D g on 

Relevant Dimensionless Parameters 

We now perform numerical computations. Our main device 
is to examine the solution of D 2 or D 2 versus a when other 
parameters are specified. Here a = 1/r) is introduced merely 
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for convenience (n.b., r\ ^ 1 thus 1 ^ a < +oo). Recall 
that when (3 ^ 1/4, is always smaller than and when 
P > 1/4, Dg can be smaller than D 2 S when a is sufficiently 
close to 1; for this reason, we show the dispersion relation 
in term of D\ through equation 1221 1 instead of Dg in most 
cases for P ^ 1/4. From now on, we use y instead of y 
to denote the D\ solution for convenience. When (3 > 1/4, 
we will show both D$ and Dg solutions with clear labels, 
respectively. 

Typically, equations 167H and 171H are algebraic cubic 
equations and have three roots and there exists at least one 
real root because all coefficients are real. In most cases, the 
three roots are all real and do not intersect with each other 
and we use yi, yi and 1/3 to denote the upper, middle and 
lower branches of the D% solution, respectively. Generally, 
they corres pond to d i fferent types of s tationary MHD den- 
sity waves (Lou 2002; Shen et al. 2005). As we have several 
parameters to choose freely in our model, the solution be- 
haviours are much more richer than previous works as we 
will discuss presently. 

As has mentioned before, the net magnetic force serves 
as either a centripetal or a centrifugal force when /3 < 1/4 
or P > 1/4, respectively. We take P = -0.24, /3 = -0.1, 
P — 1/4 and P = 0.4 in order and examine the dynami- 
cal influence of the magnetic field on solution behaviours as 
shown in Fig. 0-Fig. [I] For the moment, we focus mainly 
on the full (i.e., F = 1) disc system with m = 2. Meanwhile, 
we specify different values of mass density ratio 8 — 1/4, 
8 — 1 and 8 — 4 to assess the influence of the surface mass 
density ratio on the two coupled discs. 

The 1/1 and yi solution branches generally correspond 
to the unmagnetize d solutions [i.e., Lou & Fan 1998b; 
IShen fc Loul i2004bf) j while the 1/3 solution branch is addi- 
tional (i.e., SMDWs by nature) due to the very presence of 
a coplanar magnetic field; these MHD density wave mode 
classifications are qualitatively si milar to the correspond- 
ing isothermal solution branches of iLou fc Zoul i2004f) . This 
identification of SMDWs for j/3 branch can be seen clearly as 
the magnetic field becomes sufficiently small (e.g., q = 0.3). 
Here, the 1/3 branch is almost independent of <5 variation 
(see Fig. 0-Fig. |HJ|. In fact, this root is exactly the newly 
emerged root as noted after equation 16911 . The correspond- 
ing Dg solution of this branch remains almost constant and 
borders on zero because coefficient C$ is sufficiently small 
when magnetic field is weak and this solution will become 
trivial for q = because of zero Dg solution. 

Generally speaking, the yi branch tends to be positive 
(this may not be always true as we shall show later), and 
the 1/2 branch may be either positive or negative. The j/3 
branch is mostly negative except for few cases. The critical 
value of a when we have zero root is determined by C5 4 = 
in equation l|72|l . By a proper substitution, we find that the 
result is a cubic equation of 77 (or equivalently of a) which 
indicates that theoretically there might be three critical val- 
ues of a labelled as a c \ , &ci and <r C 3 for each branch of the 
zero root point. Practically, we often see o& and sometimes 
(j C 3 and seldom o c \ in the range of 1 ^ a < +00. However, 
for Dg solutions, since C$ is a linear function of rj, there can 
be at most only one critical value of r\ or a within a proper 
range. 

We also note that there can be situations that only 
one real root exists and the other two roots are complex 



conjugates [e.g., Fig. |l(b)| and Fig. |l(c)|. In prev i ous in - 
vestigations [e.g., |L ou fc Shen (2003~Lo u fc Zoul J2004I) : 
IShen fc Loul (|2004rf) : IShen et al.l l)2005l) ] , we never encounter 
complex roots. According to Fig. (i.e., P — —0.24), we see 
a stronger magnetic field tends to terminate the existence of 
the physical part of the yi branch. For a larger value of P, 
this tendency appears to diminish. 

As expected, the magnetic field has a very slight influ- 
ence on the composite disc system when it becomes weak 
enough (e.g., q — 0.3). For example, we compare our pane l 
(b) of Fig. El with panel (b) of Fig. 3 in lShen fc Loul <2004bl) 
and see that the y\ and yi branches are strikingly similar 
to their mutual counterparts. A point in common is that 
when a becomes large (i.e., small 77 corresponding to a sit- 
uation where the velocity dispersion in the stellar disc is 
much higher than the sound speed in the magnetized gas 
disc), the curves tend to be more or less flat. As q be- 
comes larger, or equivalently, the Alfven speed is comparable 
to or even higher than the gas sound speed, the magnetic 
field effect can be more prominent, especially for small a. 
As P becomes very small (e.g., P = —0.24), the magnetic 
field provides an attraction force, a stronger magnetic field 
should thus enhance the angular velocity of the magnetized 
gas disc (i.e., a larger Dg). Meanwhile, we see from Fig. Q 
that it reduces the angular velocity for the stellar disc. Since 
here D g > D s , so the existence of a magnetic field enlarges 
the angular speed difference between the gaseous and stellar 
discs [see also equation 1221 1. As P becomes slightly larger 
(e.g., P = —0.1), D g still increases when q grows with a 
greater attractive Lorentz force as expected (not shown in 
figures); for D s solution, the situation is more complicated. 
Here in the P = —0.1 cases with a small a, the increase of 
magnetic field strength first lowers the 1/1 and yi branches 
and then raises them upward when q becomes sufficiently 
large. One interesting phenomenon is that for large q, the y\ 
branch is flat for large a and the curve suddenly rises when a 
reduces to a certain value; meanwhile, the yi branch catches 
up and stretches that flat curve at small a [see Fig. |2(c)] Fig. 
|3(c]1 Fig. |4fc)] . For large P (e.g., P = 1/4 and P = 0.4), a 
larger q directly results in an increase in all three D 2 a solution 
branches. 

We now examine the effect of varying the surface mass 
density ratio 8. We observe a higher y\ branch for a larger 
8 = Eq/Eq and a higher yi branch for a smaller 8 when P 
and q are small [see Fig. |l(aj| and Fig. |2(a)|. This is consis- 
tent with the hydrodynamic results of lSIien fc Loul (|2004b) 
where magnetic field is absent. In the presence of magnetic 
field, we notice a new feature that as P and q becomes suf- 
ficiently large, there emerges a 'convergent point' where for 
different 8 values, some branches converge and then their 
relative positions are changed [see Fig. |3(b)| and Fig. |3(c)"| . 
For different branches, the convergent point varies. In addi- 
tion, we see from equation 170H that since X is independent 
of 8, the convergent point for Dg solution is the same as that 
for the D 2 solution. 

It is of some interest to consider the case in Fig. 0] for 
P — 0.4. Different from other cases, here D g may become 
smaller than D s in the presence of magnetic field. This hap- 
pens when a approaches 1 and q becomes fairly large. Specif- 
ically in Fig. |4(c)| Dg solution becomes negative and thus 
unphysical even though Dj. remains positive. 

Regarding the dependence on m, there are not many 
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(c) m = 2, /3 = -0.24, g = 3, F = 1 

Figure 1. Aligned solution branches derived from equations 
1671 and 1221 as functions of a = l/rj by specifying m = 2 (i.e., 
bar-like), F = 1 (i.e., full discs) and /3 = —0.24. For each set of 
<5 = l/4 (solid curve), 1 (dotted curve) and 4 (dashed curve), we 
choose q = 0.3, 1 and 3 and show corresponding results in panels 
(a), (b) and (c), respectively. In panel (a), the three curves of 
the y:j branch overlap with each other. Sometimes cubic equation 
1671 has only one real root. In panel (b), the y2 and j/3 branches 
are complex for the three chosen 5 values. In panel (c), the j/2 
and branches are complex for the <5 = 4 case For < 1/4, Dg 
remains always greater than D^. 



(c) m = 2, P = -0.1, q = 3, F = 1 

Figure 2. The solution curves of versus a for aligned m = 2 
(bar-like) full discs with F = 1, (3 = —0.1, 8 = 1/4 (solid curve), 
1 (dotted curve) and 4 (dashed curve) and q = 0.3, 1 and 2.5 in 
panels (a), (b) and (c), respectively. In panel (a), curves for the 
three different S values of the 1/3 branch almost overlap with each 
other. The same thing happens in panel (c) for large a > 1.9 in 
the j/2 branch and for small a in the y\ branch; the dashed curve 
of 3/2 branch almost coincides with the solid curve of y\ branch 
for small a < 1.9. In panel (b) with q = 1, the j/2 and j/3 branches 
are close to each other. The three curves of the 1/3 branch are 
lower than the three curves of the 3/2 branch. 
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Figure 3. The solution curves of versus a for aligned m = 2 
(bar-like) full discs with F = 1, (3 = 1/4, 5 = 1/4 (solid curve), 1 
(dotted curve) and 4 (dashed curve) and q = 0.3, 1 and 3 shown 
in panels (a), (b), and (c), respectively. Note that for small q, 
there are three physical D j. solutions as a approaches 1 and when 
q is large enough there exists a 'convergent point' where branches 
with different <5 values converge. Curves for different 5 values in 
the 2/3 branch almost merge with each other to various extents; 
in panel (a), the three curves of the j/3 branch overlap each other. 



Figure 4. The and Dg solution curves versus a for aligned 
m = 2 (bar-like) full discs with /3 = 0.4, <5 = 1, F = 1 and q = 
0.3, 1 and 2.5 shown in panels (a), (b) and (c), respectively. £>J 
solutions are shown in solid curves, while Dg solutions arc shown 
in dashed curves. The corresponding j/i, J/2i and 3/3 branches are 
ordered from up to down, respectively. In panel (c), the lowest 
Dg branch becomes negative while the corresponding D^ branch 
remains positive as a becomes sufficiently small. 
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Figure 5. The .D 2 solution curves versus a 


for aligned m 



full discs with f3 = —0.24, 5 = 1, F = 1, and g = and 1 shown 
in panels (a) and (b), respectively. For q = 0, the lower two F> 2 
solution branches actually intersect with each other. For q = 1 
and 2.3 < a < 3.2, there is no physical D 2 solutions. 



notable features except for m — 1. In this case, it is possible 
that the solution branches intersect and also there can be 
no physical D 2 solutions. We show a case in Fig. |^|with /3 = 
—0.24. We note that in panel (b) of Fig. |S] unlike previous 
m — 2 cases, here the 2/2 branch is the newly emerged branch 
due to the magnetic field. It intersects with 2/3 branch at 
first [i.e., q = in Fig. |5(a)"| . With the increase of magnetic 
field strength, it moves up and separates from the 2/3 branch 
and then meets the yi branch. The result is that complex 
conjugate roots appear and the only real root lies along the 
2/3 branch which is negative and thus unphysical. Therefore, 
the parameter regime between the two meeting points of yi 
and 2/2 branches (approximately 2.3 < a < 3.2) allows no 
stationary global MHD perturbation modes. 

3.3.3 Phase Relationships among Perturbation Variables 

Phase relationships among aligned perturbation variables 
reflect their connections on large scales and provide mod- 



els in interpreting the results from optical and synchrotron 
radio observations of barred or barred-spiral galaxies. In the 
aligned cases, most variable quantities are real (except for 
the r-component magnetic perturbation R) and thus they 
are qualitatively either 'in phase' or 'out of phase'. Our main 
interest is to examine the phase relationships between sur- 
face mass density for the two discs and their relation to the 
magnetic field. To derive a quantitative expression, we start 
from equations 1591 — 161H and obtain 



M A /(2nGE 3 Q rVr, 

S3 



@A 1 + 



(73) 



S 9 /S s = -1 
ii? mBo 

& = Hf 

Z = (4/3 - l)LR/(2m) . 

There is another version of expression for iR/S 9 obtained 
by a combination of equations l|56^ and l|57|l . namely 
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(74) 



Equations 1731 and 17411 for iR/S 3 are exactly the same for 
one can check that their combination yields the dispersion 
relation 1651 . From equation 174H . it would be easier to check 
this consistency. Setting /3 = for t he isothermal case , equa- 
tion ifTIjl reduces to equation (68) o flLou fc Zoul l|2004h (they 
considered the ratio S 9 /Z) as expected. Next is to evalu- 
ate these equations to express them in physical parameters. 
Note that iR/S 9 is no longer dimensionless, to evaluate this 
relationship, we take the proportional factor contained in 
the brackets (i.e., removing _Bo/Sq) of equation 173H for cal- 
culations. We show these results below. 
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The fir st expression i n equa tion l|75^ is the same as equation 
(42) of lShen fc Loul <2004bl) where a composite disc system 
of gravitationally coupled unmagnetized discs is analyzed. 

In equations I73H and 1741 . there is an imaginary unit 
i on R indicating that the radial component of the mag- 
netic field perturbation lags a phase 7r/2 to the mass den- 
sity perturbation. In the third expression of equation I73H . 
the factor 4/3 — 1 means that the azimuthal component of 
the magnetic field perturbation is ahead of or lags behind 
the radial component of the magnetic field perturbation by 
7r/2 for /3 > 1/4 or (3 < 1/4, respectively. In the special 
case of /3 = 1/4, we have Z — and there is no azimuthal 
component of magnetic field perturbation. 

The phase relationships for surface mass densities of 
aligned perturbations are displayed in Fig. H3-Fig- El in 
these three figures, we adopt the same parameters as used in 
computing the dispersion relations shown in Fig. 0-Fig. 21 
for complete information. For the convenience of statement, 
we still refer to the phase relation curves as yi, 1/2, and 2/3 
branches corresponding to the 2/1 , 2/2 and 2/3 of D 2 solutions 
in the dispersion relation. 

The magnetic field influences the perturbation mass 
density ratio. We demonstrate a few examples in Fig. |7| 
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Figure 6. Phase relation curves as represented by the ratio 
S 9 /S s for surface mass density perturbations versus a in the 
m = 2 cases (bar-like) for full discs with F = 1. Taking p = —0.24, 
we show curves with 8 = —1/4, 1, 4 and q = 0.3, 1, 3 in panels 
(a), (b), (c), respectively. These parameters are the same as those 
in Fig. Qfor the convenience of reference. In panel (a), the upper, 
middle, and lower curves with corresponding linetypes are the j/3, 
2/2, and y\ branches, respectively. In panel (b), only the y\ branch 
is valid and physical. In panel (c), the curves from up to down 



Figure 7. Curves of S 3 /S s ratio versus a for m = 2 full discs 
(bar-like) with F = 1, P = 1/4, 5 = 1/4, 1, 4 and q = 0.3, 3 
shown in panels (a) and (b), respectively. We adopt the same pa- 
rameters as used in Fig. |3]for reference and comparison. In panel 
(a), the curves with different linetypes from up to down corre- 
spond to 2/3, 2/2 > Vi branches, respectively, and for convenience, 
we just label the corresponding solid curves. In panel (b), there 
appear 'convergent regions' for a larger q. 



and Fig. |H] In the case of a small q (e.g., q = 0.3), the 
yi branch is negative (i.e., out of phase) while yi and 1/3 
branches are positive (i.e., in phase). However, as the mag- 
netic field increases, the y% branch becomes negative for 
small a. This means that for aligned stationary MHD den- 
sity waves, the density phase relationship between the two 
scale-free discs are out of phase for super fast MHD den- 
sity waves and in phase for fast MHD density waves. For 
the middle one, the mass density ratios are in phase in the 
regime of a weaker magnetic field but are out of phase in 
the regime of a stronger magnetic field. 

We explore phase relationships of surface mass density 
perturbations with several different 8 values in Fig. |S| and 
Fig. Q One might expect that large S leads to higher S 9 /S s 
ratio; however, this is not always the case, especially when 
both j3 and q are large (e.g., Fig. [7] and Fig. |HJ. We notice 



Coupled Magnetized Scale-Free Discs 15 



CO 

oo 





q=0.3 


/y 3 branth i . 


q=1 




q=2.5 


/ ,' m=2,p=0.4,5=1,F=1 




/ y 2 branch 




/^r-. y 1 branch ^ , ' ' 



Figure 8. Curves of S 9 /S s ratio versus a for m = 2 full discs 
with F =1, P = 0.4, <5 = 1 and <j = 0.3, 1, 2.5, respectively. The 
curves from up to down correspond to j/3, 3/2, and j/i branches, 
respectively, and for convenience, we just label the solid curve 
(i.e., g = 0.3). 



'convergent regions' where the phase relation curves with 
different S ratios seem to converge. By carefully examining 
these phase curves, they do not actually converge but just 
become very close to each other. This phenomenon results 
from a complex interaction between the magnetic field and 
the two- fluid disc system and is a new feature in our mode l 
as compared with previous works [e.g. lLou fc Shenl J2003t) : 
lLou fc Zoul <2004 : IShen fc Loul (l2004bl) ; Lou & Zou 2006]. 

We now examine phase relationships between gas mass 
density perturbation and the radial component of magnetic 
field perturbation (Fig.|5]and Fig. HUB with m — 2 (bar- like). 
From equation I74II . we know that S 3 jZ is proportional to 
S 9 /(iR) with a factor 2m/(4/3-l). When ^ 1/4 the phase 
curves regarding S 9 jZ and S 9 /(iR) are the same in shape. 
For > 1/4, S 9 /Z and S 9 /(iR) are of the same sign, while 
for /3 < 1/4, S 9 /Z and S 9 /(iR) are of the opposite sign. For 
P — 1/4, we simply have Z = 0. For these reasons we just 
discuss S 9 /(iR) for convenience; the information of S 9 /Z 
can be readily derived accordingly. In the case of p — —0.24, 
we already know that the 2/1 branch is always physical. Here 
in Fig. |3] we see that iR and S 9 remain always in phase, and 
therefore from equation 1731k we know that Z and S 9 are 
out of phase. The enhancement of the magnetic field lowers 
the ratio S 9 /(iR), as a stronger magnetic field tends to in- 
duce a larger magnetic perturbation. In the case of larger P, 
the dependence of this ratio to the magnetic field is somehow 
more sensitive. We take the case of p = 1/4 as an example of 
illustration. In this case, the azimuthal magnetic perturba- 
tion vanishes, with a nonzero r— component magnetic field 
perturbation. We see in Fig. 1101 that when q is small, the 
j/i and 2/2 branches are positive and are very close to each 
other. The 2/3 branch is negative and borders on zero (recall 
that 2/2 and 2/3 branches can be physical when a < a C 2 or 
& < o"c3). As q increases, all three branches lower and the 2/2 
branch intersects zero indicating a phase relation reversal 
across the critical point. When q is increased further, the 
entire 2/2 branch becomes negative. 
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Figure 9. Perturbation phase relation curves for gas surface mass 
density and radial component magnetic field in the m = 2 (bar- 
like) cases for full discs. Taking /3 = —0.24, we show curves with 
8 = 1/4, 1, 4 and q = 0.3, 1, 3 in panels (a), (b), and (c), 
respectively. From relation (73), we see that S 3 /(iR) and S 9 /Z 
are proportional to each other with a sign change depending on 
whether /3 < 1/4 or not. 
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Figure 10. Curves of S 9 /(iR) ratio versus it for m = 2 full discs 
with F = 1, /3 = 1/4, <5 = 1 and g = 0.3, 1, 2. The upper, middle 
and lower branches with different linetypes correspond to yi, y2 
and J/3 branches, respectively. For the solid curves with q = 0.3, 
the upper two branches y\ and yi are quite close to each other 
and the negative yz branch approaches zero. 



3.3.4 Partial Discs Embedded in a Dark Matter Halo 

Coplanar aligned MHD perturbation structures in full scale- 
free discs are extensively discussed above. For partial discs 
with factor F < 1, we expect that the introduction of an 
axisymmetric dark matter halo may lead to some novel fea- 
tures. In fact, through numerical computations we see that 
new features manifest mainly through the m = 1 case as 
shown in Fig. 1111 

We first note that one of the D$ or Dg solutions diverges 
for the cubic coefficient C % = 0. As o u r C% is proportional 
to C2 of equation (48) in IShen et alJ <l2005l) . we expect to 
have a similar result in this aspect. This phenomenon ap- 
pears at m = 1. Let /3 = 1/4 and we see in Fig. |ll(a)| the 
divergent point lies at F c ~ 0.6842 which is exactly the same 
as that in the analysis of lShen et alJ <l2005t) . and we have an 
additional 3/3 branch (although negative and thus unphys- 
ical). Another difference is that here the 3/2 branch is also 
negative and therefore the only physical solution lies in the 
positive portion of 3/1 solutions. The direct consequence is 
that for < F ^ F c , there can be no physical solutions. 

We also evaluate the m = 1 cases with different f3 val- 
ues. For a sufficiently small /3 (e.g., (5 — —0.24), there ap- 
pears a kind of discontinuity seen in Fig. |ll(b)| Actually, the 
range of F is divided into three regimes. For < F < F c i, 
only yi branch exists and the other two D$ roots are complex 
conjugates. For F c i ^ F ^ F C 2, all three branches are real. 
For F C 2 < F ^ 1, only the 3/3 branch exists. Since 3/3 branch 
is negative, there is no physical solution for F C 2 < F ^ 1. 
Immediately at F — F C 2, there emerges suddenly a new so- 
lution (double roots), and this solution then bifurcates into 
two branches as F becomes smaller. The 3/2 branch turns 
negative before F encounters F c i and from then on there is 
only one physical solution for smaller F. 

For m = 2 bar-like configurations, we take (3 = 1/4, 
8 = 1, a = 2 with q = 0.3 and 1.5 respectively to obtain so- 
lutions in panel (c) of Fig. 1111 For this set of parameters, the 
yi branch remains always positive while the 7/3 branch re- 
mains always negative. Whether the 1/2 branch is positive or 
not depends on the choice of parameters F and q. It is clear 



y,b 

- m=1,p=1/4,5=1,a=2 

F =0.6842 

\ 


ranch\ \q=1.5 

\q=0.3- N 


y 2 branch 




~~^Ty^>rahs^^ 


y tiranch 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

F 



(a) m = 1, /3 = 1/4, 5 = 1, cr = 2 



0,25 
0.2 
0,15 


"~~~~~-L m=1 
y fl branch 


,P=-0.24,5=1 ,o=1 .2,q=0.4 - 


0.1 


y 2 branch^' 




0.05 


F =0.2164 


F =0.4619 









-0.05 






-0.1 






-0,15 




y 3 branch 


-0.2 






-0,25 






0.1 0.2 0.3 0.4 


0.5 0.6 0.7 0.8 0.9 

F 


(b) 


m = 1, P = -0.24, 8 


= 1, a = 1.2, q = 0.4 


1.4 












1.2 




y 1 branch ^ 




« q-1-5 


0.8 


y 2 branch 


q=0.3 


0.6 






0.4 


m=2,p=1/4,5=1,o=2 


0.2 


y 2 branch 






-0.2 




-0.4 


y 3 branch^ 


-0.6 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



F 

(c) m = 2, /3 = 1/4, 5 = 1, cr = 2 

Figure 11. Dispersion relationships of versus F variation 
for a composite system of partial scale-free discs. There are a 
few interesting features in the case m = 1. There is only one 
real root in panel (b), while the other two roots are complex 
conjugates. It happens that the three branches are connected 
smoothly. In panel (c) for m = 2, 3/1 branch represents the upper 
two curves (q = 1.5 for dashed curve and q = 0.3 for solid curve). 
The y2 branch is labeled separately and the 3/3 branch represents 
the lower two curves. 
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that an increase of magnetic field raises all three branches. 
This raises the y\ branch only slightly but greatly lifts the 
lower two branches. Especially for 7/2 branch, only smaller F 
values can ensure a positive t/2 when q is sufficiently small, 
while a sufficiently large q parameter actually makes 1/2 pos- 
itive for all < F ^ f. For other values of m, the global 
aligned MHD perturbation patterns and trends of variation 
are qualitatively similar. 



Cj[2m 2 -if(l -4(3) + 2f 2 
2mQ. g r 2 



(f - 2\£)nGN m - 



2rE* 



+ (1 - 2iO^GAfmS s 



iU 9 = 



(1 - 2if)ag + (1 - 4j3)CT 



2Q.0 



(79) 



4 UNALIGNED MHD CONFIGURATIONS 
OF GLOBAL LOGARITHMIC SPIRALS 

In this section, we study the unaligned MHD perturbations 
by introducing a special yet useful potential-density pair re- 
ferred to as the logarithmic spirals. We first derive the MHD 
dispersion relation for the composite system of two coupled 
scale- free discs. We then focus on computations of m ^ 1 
cases. Finally, we study the m — case which leads to the 
discussion of the axisymmetric marginal stability. 



4.1 Stationary MHD Dispersion Relation 
for Global Logarithmic Spirals 

To construct global logarithmic spiral MHD perturbations, 
we take the potential-density pair in the form of 



S^r) = eV 3//2 exp(if lnr) , 
V\r) = -27rGhSWm(0 , 



(76) 



where e l is a small constant amplitude coefficient, super- 
script i denotes either g or s for either gas or stellar discs, f 
represents the radial wavenumber and the Kalnajs function 
A/" m (f) is defined by 



A/" m (f) 



F(m/2 + if/2 + l/4)r(m/2 - if/2 + 1/4) 



(77) 



2r(m/2 + if/2 + 3/4)r(m/2 - if/2 + 3/4) 

fealnaisil97fD . In potential-density pair 17611 . the perturbed 
surface mass density is complex, with an amplitude radial 
scaling of r~ 3 ^ 2 and a phase factor r' 5 for radial oscillations. 
The relevant radial variations are no longer the same as 
those of the background equilibrium unless (3 = 1/4. In fact, 
th ere can be m o re gene ral potential-density pairs as noted 
in IShen fc Loul i2004bT) . For the analysis of the marginal 
stability problem, expression I76H is used such that the dis- 
persion relation 198H derived later is real on both sides, and 
this choice of scale factor satisfies the scale-free condition in 
the perturbation equations. The phase factor is the key to 
establish global logarithmic spiral patterns. We call it loga- 
rithmic spirals because the dependence of the phase factor 
on r follows the logarithmic function. Kalnajs function 17711 
is an even function in f and thus a consideration of f ^ 
(i.e., leading MHD density waves) suffices. 

We now derive the MHD dispersion relation for global 
logarithmic spirals with m ^ 1 cases using equations 14211 — 

dlH . together with With radial scalings of 5" oc 

r -s/JH-tt v x r -i/a+« u x r -i/2+/3+iC and j x r i/2+^+i 5j 

all quantities in equations 142H — I44H and 1481 1 are scale- free; 
these equations become 



mn g S 9 - (1/2 - if + (3)%U 9 + ^J 9 = , 



(78) 



(l-/9)rn 9 



2nGN„. 



Ci(l-4/3) 



2rfi„ 



iU 9 



rEg 



S 9 + 2nGmrAf m S 3 



Z = -i£ 



B iU 9 



R = 



BpU 9 



for the magnetized gas disc and 



(m 2 -2 + 2(3)Q? 3 r 



+ f 2 + 2/3 + l/4 



- a 2 + 27rGrEo7V„, 



■ 2nGrY, s Af m S 9 = 



(80) 



(81) 



(82) 



for the stellar di sc. Equatio ns 178[ — Q81I are the same as 
equation (65) in lShen et alJ ([2005) for a single magnetized 
scale- free disc by setting S s = and equation 1821 is ex- 
actly th e same as the first equation of (68) in lShen fc Loul 
<2004tft . A combination of equations CHI - JHU and d%2l 
gives the full solution for the stationary dispersion relation 
of the gravitational coupled MHD discs for logarithmic spi- 
rals. Similar to the aligned cases, these are homogeneous lin- 
ear equations in terms of (S 9 , S s , U 9 , J 9 ). For non-trivial 
solutions, the determinant of the coefficients of equations 
178t — 1801 and 18211 should vanish. The result is the station- 
ary MHD dispersion relation. As necessary requirements, we 
expect that in special cas es, this dispersion rel ation will re- 
duce to those discussed bv lShen et alJ ll2005l) or lShen fc Loul 
(2004b). Meanwhile, we solve the above equations directly 
to reveal physical contents. 

A combination of equations I78K and 1801 leads to two 
expressions of iU 9 and J 9 in terms of S 9 and S s , namely 



iU 9 = mQ g r[(@ s + E S )S 9 + 6 S 5 S ]/Eg , 

J 9 = 0,r 2 {[-l + (1/2 - if + 0)(Qs + E S )]S 9 
+ (1/2 -if + /3)9 s S a }/Eg , 

where, for notational simplicity, we denote 

6s = 2nGrAf m Y, 9 /A s , 
E S = (-a 2 + fi 2 r 2 )/A s , 
A s = (3/2 - if )Q 2 g r 2 + (2(3 - 



(83) 



(84) 



l/2)Cl 



We use the subscript s to indicate the logarithmic 'spiral 
cases'. Note that Os and Es are two dimensionless constant 
parameters. A substitution of equation I83K into equation 
1791 1 and a further combination with equation 1821 1 lead to 



[Ks(6s + E s ) +jC s - (1 - 2iti)nGY, 9 rAr m ]S 9 

+ [/C s 6s - (1 - 2if)7rGEgrP m ]5' s = 
2TTGY, a rAf m S 9 + (Ms + 2nGZ s rAf m )S s = 



(85) 
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where, for notational simplicity, we further define 

Ks = (m 2 + 2/3 + 1 - 2i^)Q 2 s r 2 

-[m 2 -i£(l/2-2/3)+£ 2 ]Ci , 
C s = (1/2 - 2f3)C 2 A - 2Q, 2 g r 2 + (1/2 - U-)a 2 g , ( 86 ) 
(m 2 + 2/3 - 2)Q, 2 r 2 2 



The four coefficients Cf , C* 2 S , C*f, C| are functions of di- 



Mi 



(m 2 + e 2 + 2/3 + 1/4) 



We can derive the stationary MHD dispersion relation for 
global logarithmic spirals by calculating the coefficient de- 
terminant of equations l|82|l and l|85[l . After manipulations 
and rearrangements, we obtain 

{A s /Ks)[lCs(@s + S s ) +£ g -(l- 2iC)7rGEgrA^ m ] 

x {Ms + 27rGEgrA/ m ) (87) 
= 4^ 2 G 2 E«E^ 2 AT 2 [1 - (1/2 - i0As/fcs] . 

In equation 1871 . the left-hand side consists of two factors. 
One can show that the factor in the left brackets is ex- 
actly the dispersi on relation in a sing le coplanar magnetized 
disc discussed by IShen et all i2005lf and the second factor 
in parentheses denotes the dispersion relation of a single 
scale-free stellar disc. The right-hand side denotes the mu- 
tual gravitational coupling between the two scale-free discs. 
By setting j3 = for the i sothermal case, eq uation I87H re- 
duces to equation (106) in lLou fc Zoul i2004T) where isother- 
mal fluid-magnetofiuid discs are investigated. 

In all calculations of this subsection, the equations for 
logarithmic spirals are strikingly similar to those of the 
aligned perturbations. If we replace V m by M m and partly 
replace 2/3 by 1/2 — i^, this similarity arises in the expression 
of the perturbation equations. We assert that the similarity 
in structure must lead to the similarity in concrete expres- 
sions of the dispersion relation between the two cases. We 
show this fact in the next subsection. 



4.2 Stationary Logarithmic Spiral Configurations 

4-2.1 D 2 and D 2 Solutions of the Dispersion Relation 

We now evaluate the spiral dispersion relation for m > 1 
cases numerically to show the dependence of D 2 g or D 2 on 
other dimensionless parameters (i.e., m, /3, 77, 5, q, F). 
Similar to the aligned cases, we choose to calculate the co- 
efficient determinant of equations 178H — (1801) and 1821 . We 
leave details of onerous calculations to Appendix C where 
the simplified determinant is shown explicitly. After careful 
calculations and simplifications, we obtain the MHD disper- 
sion relation as cubic polynomial equation in terms of D 2 . 
It is satisfying to see that this dispersion relation can be 
expressed in the exactly same form as that of the aligned 
cases merely by a substitution of notations. We write 

aUp, = m 2 + <e 2 + 2/3 + 1/4 , 

H s m (P, 0=A s m M m /C + B m . 

We still use the notations B m (P) and C(/3) as defined in 
equation 16611 . Now the stationary MHD dispersion relation 
can be written in terms of y = D 2 namely 



Ciy 3 + Cly 2 + Cfy + C$ 



mensionless parameters m, /3, 5, 77, and q 2 defined by 
{Af n + B m ){H s m -B m )b- 



c! 



(1 + 5) 



(Ai + B m )(H? n S + B m ) s 2 
r ^2QQ 



A s T-/ s 4- 



(1 + 5) 
{A s m + B m ){H s m 



B m )5 



(1 + 5) 
(Ai + B m )(Hi5 + B„ 



(90) 



(1 + 5) 

iS 2 



Cf Q q 2 



+ (A^ + Bm) 

Co = C0Q3? 6 + C*oq29 4 + CoQiq 2 
where the following five coefficients are further defined by 
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Ciq 
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2AZ, + I 
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4(1-2/3)^ - (1 + 2/3) B m M 



20/3 + 6 , 

{rim + BmOj-q 



(4/3 + 2)(l + 5) 
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■Sri , 



A s m + 5/3 - 



(1 + 5) 

^(4/3-3) + 2/3+l/ s 2B m N m 8 \ 

TL m + DmO Tj 



(4/3 + 2)(l + 5) 
_ ^lg t (4/3-3) + 2/3+l A' Tr , 
(4/3 + 2)(l + 5) C 

+ (Ai + B m )(Ai + 5(3 - 2) 



{A b m + B m )5 



^0Q3 



+ 



(4/3 -l) 3 (H~l + BmS)y 

2(4/3 + 2) 2 (l + 5) 
(4/3-l) 2 (2^-2/3-l)g m AC 

C 

2/3- 



(4/3 + 2)3(1+5) 
(4/3-1 



~ (4/3 + 2) 

x {A s m + B, 
(4/3 - l)(2Ai 



(4/3 - 1) 2Al 



Si] , 
- 1 M m 



(4/3 + 2)2(1 + 5) 



2 ' (4/3 + 2)(l + 5) C 

(4/3-l) 2 (^g t + g m 5)? ? 
2(4/3 + 2)(l + 5) 

2/3 - 1) f„s j_ n x 2B m M m S 
rt m + o m o -= 



<^0Q1 



(2A s m - 2/3 - 1) 
= (4/3 + 2)(l + 5) 

(2A s m - 2/3 - 1) 
(4/3 + 2)(l+5) 



Tit + Bm5 



B m M m 8 
C 



l + 5-^S)(Ai--B ll 



(91) 



(89) 



Before solving dispersion relation l|89jl . we see qualitatively 
the dependence of y = D 2 on other dimensionless param- 
eters. Similar to the aligned cases, we have Ci = for 
q — and meanwhile, C§, Cf, Cf are exac tly the same 
as C2, Ci, Co respectively in equation (72) of IShen fc Loul 
(2004bj). By setting r\ — and 5 — > 00, we have C3 = 
and the coefficients actually reduce to the situation of a 
single MHD scale-free disc analyzed by IShen et al.1 (120051) 
[see their equations (70)]. We also point out that by setting 
/3 = 1/4 and £ = 0, we have A/" m (l/4, 0) = P m (l/4) and 
-4m(l/4, 0) = .4^(0), this logarithmic spiral case goes back 
to the aligned /3 = 1/4 case as expected. 
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Equation 189H is cubic in y = Dg. We can also readily 
cast equation (I89H in terms of Dl instead of Dg. This is 
exactly parallel to what we have done in the analysis of the 
aligned cases and we would not repeat the steps here. For 
notational clarity, we should add ' to the coefficients of the 
cubic equation for D 2 S such as Cq and so forth, similar to 
the aligned cases. However, since Dl will be frequently used, 
we denote y = Dl instead of y' from now on unless otherwise 
stated. 

4-2.2 Dependence of Dl Solutions with 
Other Dimensionless Parameters 

As we have shown before, the dispersion relation for global 
logarithmic spirals is expressed in the same form as that 
for aligned perturbations, indicating that the behaviours of 
the two cases bear certain similarities. Compared with the 
aligned cases, the differences should arise only from the sub- 
stitutions of At(P) by A^(J3, f) and V m (f3) by Af m (0, re- 
spectively. For logarithmic spiral perturbations, there is one 
more parameter £ corresponding to the radial wavenumber. 
This parameter is important for forming large-scale struc- 
tures in spiral galaxies. On the one hand, we are interested 
in the influence of £ on behaviours of Dl solution. On the 
other hand, by analogy we can imagine that the dependence 
of Dl solution on parameter r\ should be similar to that of 
the aligned cases. For these two reasons, we present our fig- 
ures in £ — y plane where y = Dl. We focus on the m = 2 
case (viz., two-armed spirals). For cases of different (3 values, 
we will show in Fig. I12t -Fig. 1141 The m — 1 case is shown 
in Fig. 1151 briefly. We also show results for partial discs in 
Fig.Gl 

Similar to the aligned cases, the yi and 1/2 branches cor- 
resp ond to the unmagnet ized solutions of Lou & Fan (1998b) 
and IShen fc Loul 1120040) and 1/3 is an extra solution (i.e., 
SMDWs) due to the very existence of magnetic field. Mean- 
while, the upper two branches are monotonically increasing 
with the increase of £, while the j/3 branch remains always 
flat as £ varies. The 2/3 branch is sensitive to the q variation. 
This is a sign that the 1/3 solution is the direct result of the 
magnetic field. Generally, the yi branch remains always pos- 
itive, the 1/2 branch can be either positive or negative and 
the lowest j/3 branch is mostly negative and thus unphysical. 
There are critical values of £ at which these branches change 
signs. These critical values are determined by Gq — 0. As 
they cannot be evaluated analytically, we simply show their 
presence in the figures. In each figure, we choose a — 1 and 
a — 5 and see that all three branches are lowered for larger 
values of a as expected. This trend of variation is similar to 
that in the aligned cases. 

For logarithmic spiral cases with m = 2, the 'convergent 
point' is a common phenomenon where curves with different 
S values tend to converge at a certain point. For a small (3 
(e.g., (3 — —0.24), there is no convergent point in the aligned 
cases but here the variation of £ produces this convergent 
point. For larger f3 values (e.g., (3 = 0, 1/4), the situation 
becomes similar to the aligned cases where variations of a 
lead to the existence of convergent point. For example, in 
Fig. ESI and Fig. HI the sequence of curves with different S 
values show a sign of convergence when a — 5. 

A strong coplanar magnetic field generally raises all 
three branches [see Fig. |13(b)| and Fig. 1141 . in accordance 
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Figure 12. D 2 solution branches of equations 1891 and 1221 

as functions of radial wave number § with specified parameters 
m = 2 and (3 = —0.24. When q = 0, our cas es naturally reduce to 
the unmagnetized coupled discs analyzed bv lShen fc LovJ J2004b) 
where there are two roots corresponding to faster and slower den- 
sity waves (Lou & Fan 1998b). For the convenience of comparison, 
we draw this case of q = in boldface solid curves. When a = 1 
and 9 = 0, the solution is independent of S JShen fc Lorjl2004bl) . 
Setting q = 2, we then show the solutions with 8 = 1/4, 1 and 
4 which are displayed in solid, dashed and dash-dotted curves, 
respectively. We examine full discs with a = 1 and a = 5 shown 
in panels (a) and (b), while for a partial disc, we show a case with 
(7=1 and F = 0.5. 
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Figure 13. D'j solution curves versus £ for a logarithmic spiral 
case with m = 2, /3 = 0, F = 1. Magnetic field has a considerable 
influence on the disc system [see panels (a) and (b)]. A weak mag- 
netic field results in a physical 2/3 solution but when the magnetic 
field strength becomes stronger, this 1/3 solution falls negative. In 
panel (c), we show the influence of cr, where the three curves of 
j/3 branch pact together. 



Figure 14. solution curves versus £ for a logarithmic spiral 
with m = 2 (two-armed spirals) and (3 = 1/4. By setting q = 1.5, 
we show the solutions with <5 = 1/4, 1 and 4, respectively. The 
cases for full discs with u = 1 and u = 5 are shown in panels (a) 
and (b), respectively. For a partial disc, we set a = 1 and F = 0.5 
and show the results in panel (c). The three curves of branch yz 
overlap with each other. 
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Figure 15. The Z)J solution curves versus £ for a logarithmic 
spiral with m = 1. For /3 = —0.1 in panel (a), there can be only 
one real root and for /3 = 1/4 in panel (b) the curves are similar 
to the unmagnetized case. In panel (b), the three curves of j/2 
branch almost overlap with each other. 



with £. For a small 0, there can be only one real root and 
the other two are complex conjugates [see Fig. |15(a)| . We 
also see that both q and S have a strong influence on the 
solution behaviours. For the physical yi solution branch, a 
strong magnetic field as well as a high 5 greatly raises the 
solution curves. 

For a composite system of two partial scale-free discs 
embedded in an axisymmetric dark matter halo, we distin- 
guish m — 1 and m > 1 cases. In the m = 1 case, there 
always exists the critical value F c for the potential ratio F 
parameter such that the j/i solution branch diverges from 
either side of F c . Mathematically, this critical value F c is 
determined by C§ = and is quite similar to that in the 
aligned cases. For F ^ F c , there is no physical solution, 
while for F > F c the yi branch is physical and sensitive 
to F variation; the existence of an axisymmetric dark mat- 
ter halo can significantly influence global configurations in 
a composite disc system with m = 1 MHD perturbations. 
Here F c depends only on /3 and £ and we demonstrate this 
dependence using examples in Fig. |16(b)| The F c curves 
with different j3 decrease monotonically as £ increases. For 
a small radial wavenumber £ ^ 1, the critical value F c ap- 
proaches 1 as seen in panels (a) and (b) of Fig. 1161 Since 
in reality the dark matter halos play significant roles in disc 
galaxies (i.e., F is small and not close to 1), we propose that 
a spiral galaxy with only one spiral arm with a relatively 
large proportion of dark matter halo must have a relatively 
large radial wavenumber (i.e., tightly wound) to support a 
global stationary pattern. In principle, this may be tested 
by observations. 

For spiral perturbations with m > 1, however, the dark 
matter seems to be less important in determining the sta- 
tionary dispersion relation. Because in the entire range of 
< F ^ 1, there are only slight variations in the D 2 S solu- 
tions. We show such an example in Fig. |16(c)| With a radial 
wavenumber f = 2, all three solution branches j/i, 2/2 and 
j/3 remain positive and rise only slightly as F decreases. We 
can also compare Fig. |12(a)| and Fig. |12(c)| Fig. |14(a)| and 
Fig. |14(c)| and see that 'mm = 2 cases, the dark matter halo 
has little influence on solutions. For cases of m 3, the situ- 
ations are qualitatively similar. Nevertheless, a dark matter 
halo should have a strong stabilizing effect in general. 



with the aligned cases. Nevertheless, this is not always the 
case. For a small /3 (e.g., /3 = —0.24), a larger q tends to 
lower the solution curve as the radial wavenumber £ be- 
comes small (see Fig. 1121 . In the (3 — case corresponding 
a composite system of two coupled SIDs (Lou & Zou 2004, 
2006), the increase of magnetic field strength first lowers the 
j/i and 2/2 branches for small £ and then turns to raise them. 
Meanwhile, 2/3 is positive when q is not too large, but as 
q grows, 1/3 is lowered down to become negative and thus 
unphysical (see Fig. I13H . We made a few comments about 
the 2/3 branch. We already see that it bear little relation 
to the spiral pattern. Since it is the new branch originated 
from the magnetic field, we have 2/3 = when q = and its 
dependence on q is easily seen from our figures. Generally 
speaking, the 2/3 branch is sensitive to both a and q. 

The m — 1 case reveals some interesting features to 
note. We show a few cases in F ig. |l5l Here the unmagnetized 
two solutions of IShen fc Loul (l2004bf) correspond to our yi 
and 2/3 branches. The 2/3 branch is monotonically decreasing 



4-2.3 Phase Relationships among Perturbation Variables 

For logarithmic spiral cases, the phase relationships among 
coplanar MHD perturbation variables have more features to 
explore. The introduction of radial oscillations with radial 
wavenumber £ brings an imaginary part into our equations 
which describes the phase difference between perturbed sur- 
face mass density and the magnetic field perturbation as well 
as other perturbation variables. To describe these phase re- 
lations more specifically, we begin with equations 182H — 1841 
to obtain 

S 3 /S a = -1 - Ms/{2-KGY, s rN m ) , 

m mB r / s s \ i 

^=idH i+ sF)H ' (92) 

Z = £LR/m . 
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(solid curves) and 1.5 (dashed curve) 
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(b) Critical value of F c versus £ for the case m = 1 
with three different (3 values. 



There is another independent expression of iR/S 9 obtained 
by a combination of equations 1781 and 17911 . namely 



iR 

S3 



mBo 
ys 

Ml 



(1 - 2iQ7rGSgrjV m 
ICs 



1 + 



S3 



ICs 



(93) 



It is no wonder that these expressions are quite similar to the 
aligned phase relationships. We next express these results in 
physical parameters. As iR/S 9 is no longer dimensionless, 
we take the proportional factor in brackets of equation 19211 . 
The results are 



S" 

m 

S3 
+ 



(B m D'i-Ai)C(l + S) 



(£>§ 



(l + 2/3)(3/2-i^-(l/2 
D 2 g + 1 + (4/3 - l)g 2 /(4/3 + 2) 



C(l + S) 



W)i 2 
s- 



(l + 2p)D 



(94) 



Since iR/S 3 contains imaginary part, the phase difference 
between surface perturbation mass density and the per- 
turbed magnetic field becomes complicated. For this reason, 
we show such phase relations in special figures. 

We first quickly note that when q — (unmagnetized) 
and cr = l (the same 'sound speeds'), the ratio S 9 / S s = ±1 
for the j/i and j/2 solution branches, respectively. We can 
prove this claim by using the fact that yi = Am/B m 
and V2 = (1 — Mm lC)Am.l'Hm for the two D 2 solutions 
iShen fc Loull20 04b') and by substituting these explicit y\ 
and 7/2 expressions into equation 11941 . For fairly arbitrary 
parameters, the tendency of the phase relation curves for the 
ratio of the perturbed surface mass densities in the two scale- 
free discs are shown in Fig.f^and Fig.fK] We would empha- 
size that the relative order of the three branches varies in the 
phase relation curves. For example, for a small f3 = —0.24, 
the j/3 branch is the lowest, while for a larger (3 — 1/4, it 
becomes the upper one. 




(c) to = 2, p = 1/4, 5 = 1, er = 2, £ = 2, q = 0.3 
(solid curves) and 1.5 (dashed curve) 



Figure 16. Solutions for a composite system of partial scale-free 
discs. For m = 1 cases, there exists F c such that the y\ solution 
diverges [see panel (a)]. This F c depends only on P and £ as shown 
in panel (b). For m ^ 2 cases, the ratio F does not have much 
influence on the behaviour of D 2 solutions. 
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Figure 17. Phase relation curves for the ratio S 9 /S s of surface 
mass density perturbation versus the radial wave number £ for 
the case of m = 2 and J3 = —0.24. The solutions with q = 3, 8 = 
1/4, 1 and 4 are shown in light solid curve, dashed curve and 
dash-dotted curve, respectively. We examine full discs (F = 1) 
with ct = 1 and cr = 5 in panels (a) and (b), respectively. For a 
partial disc with F = 0.5, we show solutions with cr = 1 in panel 
(c). In all panels, the solution branches y±, 1/2, and ys can be 
readily identified. Each branch contains three distinct solutions 
with (5 = 1/4 (solid curve), 1 (dashed curve) and 4 (dash-dotted 
curve), respectively. In each panel, the two boldface curves for y\ 
and j/2 branches correspond to q = (the ys branch is trivial). 



(b) m = 2, f3 = 1/4, cr = 3, F = 1 



CO 





S=1,q=0 




S=1/4,q=3 


y 3 branc^i^^-^^^^^^^^ 


- - 5=1,q=3 




- - 8=4,q=3 










y ? branch 






m=2,p=1/4,a=1,F=0.5 


y fl branch 




y fl branch 







1 2 3 4 5 6 



(c) m = 2, f3 = 1/4, cr = 1, F = 0.5 

Figure 18. The phase relation curves for the ratio S 3 / S s of sur- 
face mass density perturbations versus £ for a logarithmic spiral 
case with m = 2 and (3 = 1/4. We show separately phase curves 
for full discs (F = 1) with cr = 1 and cr = 3 in panels (a) and (b), 
respectively. For a partial disc of F = 0.5 with cr = 1, we show so- 
lution branches in panel (c). Note that in the case q = (boldface 
curves) and cr = 1, t he two solution bran ches have S 9 /S s = ±1 
[see equation (73) of lShen fc Lol <2004bM . 
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In the case of /3 = —0.24 and for full discs, the phase 
relation of the perturbed surface mass densities in the two 
scale- free discs is out of phase for the yi branch and in phase 
for the 1/2 branch (yz branch is unphysical) even with a 
strong magnetic field; this situation differs from the aligned 
cases. For the yi branch in full discs, a small £ leads to 
a larger S 9 / S s ratio, while for a large £, the ratio can be 
rather small compared with 8. For the yi branch in full discs, 
the S 9 1 S s curves with different 8 values do not vary signif- 
icantly. For a small £, the ratio \S 9 /S"\ is fairly small and 
for a large £, this ratio lies between 1 and 2. The increase of 
a generally raises the two branches to different levels. The 
presence of magnetic field gives rise to variations in the y\ 
solution curves in two trends: the curves move above and 
below the q — solution curve for small and large £, respec- 
tively (see Fig. 1171 . For a composite system of two coupled 
partial scale-free discs, some of these tendencies does not 
hold. First, S 9 /S s may become positive for the yi branch 
with a sufficiently small radial wavenumber £. It also raises 
the t/2 branch considerably and separates solutions of y\ 
branch with different 8 values. In other words, the influence 
of 8 becomes more prominent. 

In the case of /3 = 1/4, the phase relation curves of the 
perturbed surface mass densities for the yi solution branch 
remain negative. For small a, variations of 8 produce con- 
siderable differences in phase curves, while for large <r, phase 
curves with various 8 values become much closer. For the 7/2 
branch with a small a, the ratio S 9 /S s can become negative 
at a certain a c with S 9 = and the curves with different 8 
values are bunched together. For a large a, the ratio S 9 /S s 
are positive and the curves with different 8 values are sepa- 
rated. In the case of a small a = 1, the yz branch is physical 
with a positive S 9 jS s ratio growing with £. The value of 
8 influences the phase relation in this branch in a signif- 
icant manner. In a composite system of partial discs, the 
presence of a dark matter halo raises the 1/2 and yz solution 
branches and lowers the y\ solution branch in the phase re- 
lation curve. Meanwhile, the S 9 jS s ratio is positive for the 
j/2 branch, and for the y\ and yz branches, the ratio magni- 
tude \S 9 /S"\ becomes much larger than that for a composite 
system of full discs. 

We now examine the phase relationship between the 
perturbed gas surface mass density and the perturbed mag- 
netic field as these aspects may provide useful diagnostics for 
observations. Since both S 9 /(iil) and S 9 /Z are complex in 
general, we show such phase relation in the complex plane. 
In Fig. E3 and Fig. [Xfi for the complex ratio S 9 /(iR), the 
variation of £ is labelled by an asterisk * with an interval of 
A£ = 0.3 starting from the real axis. According to equation 
19411 . the S B /(iR) ratio is real for £ = and thus the curves 
start from the real axis. In Fig.|5Tland Fig.|Hfor S 9 /Z, the 
variation of £ is labelled by an asterisk * with an interval of 
A£ = 0.4. According to equation S 9 /Z — > cx) as £ — > 0. 
Equation 192H also tells us that the phase relations of the 
azimuthal magnetic field perturbation and the gas surface 
mass density are totally opposite for leading (£ > 0) and 
trailing (£ < 0) spiral MHD density waves. 

In the case of /3 = —0.24, only the y\ and 1/2 branches 
are physical, and we show these two branches in Fig. ll9l and 
Fig. 1211 The main trend is that with the increase of the ra- 
dial wavenumber £, the real part of S 9 /(LR) is positive and 
does not vary much as the imaginary part decreases steadily 
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Figure 19. The phase relation curves for the perturbed radial 
magnetic field and gas surface mass density versus £ for the case 
of m = 2 and [3 = —0.24. The ratio S 9 /(i-R) is generally complex 
and is shown here in the complex plane with real and imaginary 
axes. We take the range of £ = to 6. For £ = 0, the imaginary 
part is zero according to equation 1941 . As § grows, the imaginary 
part emerges and becomes more significant. To demonstrate the 
relation with £, we show * to mark an interval of 0.3 in the vari- 
ation of £ from to 6. In this case, only the yi and 1/2 solution 
branches are physical and the unphysical yz branch is ignored. 
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Figure 20. The phase relation curves for real and imaginary 
parts of S 9 /(iR) between the radial component of the perturba- 
tion magnetic field and surface mass density versus £ for loga- 
rithmic spirals with m = 2 and = 1/4 in a complex plane. We 
show separately curves of full discs (F = 1) with cr = 1 and a = 3 
in panels (a) and (b), respectively, and for a partial disc with 
cr = 1 in panel (c). We also add the phase curve for the physical 
yz branch. When q = 1.5, the yz branch is very far away from the 
other curves; we have omitted this curve for the compactness of 
presentation. 
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Figure 21. The phase relation curves between the perturbed 
azimuthal magnetic field and surface mass density versus cr for 
logarithmic spirals with m = 2, = —0.24. We show curves of 
full discs (F = 1) with cr = 1 and cr = 3 in panels (a) and (b). 
We take the range of £ = 0.4 to 4. Since S 9 /Z — > oo as £ — > 0, 
we start from £ = 0.4 and use * to label an interval of 0.4 in 
the variation of £. The phase points moves from right to left as £ 
increases. 



from zero. The result is that for tight-winding spiral arms 
with large £, the radial component of the magnetic field 
perturbation becomes almost in phase with the gas density 
perturbation (i.e., S 9 / R — » a positive number), while for a 
small £, the radial component of the magnetic field pertur- 
bation lags behind the perturbed gas density by ~ n/2, a 
trend similar to the aligned cases. The azimuthal compo- 
nent of the magnetic field perturbation Z = £iR/m remains 
always ahead of the radial component of the magnetic field 
perturbation by 7r/2 in phase. We note that for small £ the 
imaginary part of S 9 jZ approaches a constant as its real 
part moves to infinity. Meanwhile, the azimuthal magnetic 
field perturbation can be much larger than the radial mag- 
netic field perturbation for large £• For a composite system 
of partial scale-free discs with a coplanar magnetic field, 
we note the following two points. First, by decreasing the 
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Figure 22. The phase relation curves between the perturbed 
azimuthal magnetic field and surface mass density versus a for 
logarithmic spirals with m = 2, /3 = 1/4. We show curves of full 
discs (F = 1) with a = 1 and c = 3 in panels (a) and (b). We also 
add the physical j/3 branch. When q = 1.5, the 1/3 branch is very 
far away from the other curves and we have omitted this curve 
for the compactness of presentation. We start from § = 0.4 and 
use * to label an interval of 0.4 in the variation of £. The phase 
points moves from right to left for the y\ and j/2 branches and 
from left to right for the 1/3 branch as £ increases. 



potential parameter F, the 2/2 branch moves towards right 
while the j/i branch does not move much. Secondly, for a 
large magnetic parameter q = 3, singularity occurs in the yi 
branch as F decreases. This singularity can be determined 
from equation (1941 . namely 



[D 2 + 1 + (4/3 - l)g 2 /(4/3 + 2)]Af m 5 



C(l + S) 



, s 3 

S9 



(95) 



-1 + (1 + 2/3)D 2 B 







where F parameter is implicitly contained in coefficient C. 

The 2/3 branch can be positive in the case of /3 = 1/4 
(see Fig. 1141 for this SMDW mode) ; we therefore include 
the phase information of the 1/3 branch in Fig. 1201 and Fig. 
1221 Different from the 1/1 and 3/2 branches, the 1/3 branch 



has a S 9 /(iR) characterized by an almost negative real part 
and an imaginary part increasing (rather than decreasing) 
steadily with the increase of £. The phase relation of this 1/3 
branch has a phase difference of ~ ir relative to the other 
two branches. For S 3 /Z, the 1/3 branch still lies in the sec- 
ond quadrant and moves from left to right and tends to the 
upper imaginary axis. It indicates that for this branch the 
azimuthal magnetic perturbation lags behind the gas sur- 
face mass density perturbation by ~ n for small £ and by 
~ 7r/2 as £ becomes fairly large. There is yet another salient 
feature. For a small a = 1 and a large q = 1.5, the 1/2 branch 
appears fairly special. In panel (a) of Fig. I2UI starting from 
a negative real value as £ increases, the imaginary part first 
increases and then decreases, meanwhile the real part in- 
creases with increasing £; the curve intersects the real axis 
again at the origin exactly. It then keeps this tendency such 
that the imaginary part falls down steadily as £ increases. 
In terms of phase difference between S 3 and R or S 3 and Z, 
this branch is special in that there is a phase shift of n at 
some £ c . This point corresponds to S 3 — 0, also related to 
the phase shift of S 3 / S s discussed earlier. For cases of par- 
tial discs with F — 0.5, such 'singularity' no longer exists 
and all curves are continuous. Meanwhile, such ~ ir phase 
shift of y2 branch in the presence of a strong magnetic field 
no longer exists either. 



4.3 Marginal Stability Curves for 
Axisymmetric Disturbances 

The m = 'spiral' case is very special. It is considered as 
a purely radial oscillation in a composite MHD disc system 
and is related to part of the MHD disc stability problem. 
In a general disc problem, the dispersion relation for time- 
dependent MHD perturbations involves a to 2 term. A lin- 
early stable disc system requires cu 2 ^ 0. In our analysis, we 
set lu — to obtain stationary patterns and to also meet the 
requirement of the scale-free condition. In the m = case, 
the stationary dispersion relation with u i = actually cor- 
responds to the marginal stabil i ty cur ve (IShen fc Loull2003l. 
1200(1 l2004bl : iLou fc Zoul 120041 . 12006ft . Meanwhile to treat 
the problem properly, a different limiting process should be 
taken similar to the aligned m — case. That is, we first 
set m = with uj 7^ and then let uj — > in the rel- 
evant MHD perturbation equations. The same as before, 
we begin with equations 14911 — 1531 . Equation 1491 requires 
U oc r 1 /2+2 / 3+ij and then equation (gJJ gives J oc r 1/2+f3+ ^. 

Substituting these two expressions into equation I50L we 
see that the scale-free condition cannot be met unless u; = 
where U — with U/lo 7^ 0. In the process of taking the 
limit of lu — * 0, we also remove terms that are not scale free. 
We show the result of this limiting process below. 



r 2 (iu 2 - n 2 g )S 3 = + l/4)Sg* 9 + (C 2 + 1/4 - I3)C 2 A S 3 
t 2 {lo 2 - k 2 s )S 3 = (f + 1/4)ES* S . 



(96) 
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A substitution of equations 1361 and 17611 into equation 19611 
with uj — yields 



[AS a 2 g + n 2 g r 2 + (AS - I3)C 2 A - 2-KGrAZ A/" Eg]S 9 

- 2tyGtAo VoEgS 8 = 
[Aq fls + K 2 r 2 — 2-irGrAo J\foT>o]S s 



(97) 



2-kGtAq 'Mo^S 9 







where we define 

At' =e+ 1/4 



and the corresponding notation 
Hi' = Al'Afo/C + Bo . 

The difference between Aq and A$ is that the latter does 
not contain the term 2/3; this is exactly the result of different 
limiting processes. 

The two relations in equation 1)9 7|l cannot be satisfied 
unless the coefficient determinant vanishes, which leads to 

[Al'a 2 . + n 2 a r 2 + (A§' - /3)C 2 A - 2ivGrA$' .A/bE 9 ,] 

x {A$' a 2 + K 2 r 2 - 2-kGtAI' N Y, S ] (98) 

= {2-kGvAS ' Mof^l . 

Substitutions of equations (121 )t — i'l'M lead to the explicit 
expression of this stationary dispersion relation 198H in terms 
of physical parameters in the form of 
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By rearranging equation 1991 . we obtain a quadratic equa- 
tion for the marginal stability in terms of y — D 2 , namely 
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Figure 23. The marginal stability curves for the case of m = 0, 
(3 = —1/8, 6 = 1, cr = 5. For comparison, we show the cases q = 
(boldface curves) and g = 1.5 (light curves) in one figure. Note 
that for a full disc (F = 1) and a partial disc (F = 0.6), the 
marginal stability curves behave quite differently. 
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One necessary check for a consistent derivation is to set 
q = and the coefficients <f 10 1 |l reduce to those in the case o f 
unmagnetized coupled scale- free discs (Shc n fc Loull2004rJ) . 
Equation 1100H is quadratic in D 2 and has two solution 
branches corresponding to the marginal stability curves. For 
several chosen sets of parameters, we demonstrate these D 2 
curves in Fig. 1231 and Fig. 1241 The unstable regime consists 
of two regions. In the upper-right corner, the two coupled 
discs rotate too fast to be stable and are susceptible to the 
MHD ring fragmentation instability, while in the lower-left 
corner the two coupled discs rotate too slow to resist large- 
scale Jeans collapse modified by rotation and magnetic field 
jLernos et al.lll99lt ISver fc Tremaindll99dlShu et alJl2000l: 



Fanll2002l: lLou fc Shed \200A lLou fc W 



Vu 
)4, 



Lod|2002l: lLou 

2005[lShen fcLoul2004bl:IShen et aU2005tlLou fc Zoul20Q. 

2006). The marginal curves represent precise global treat- 
ments of axisymmetric linear stability. In the WKBJ or 
tight-winding regime, one can readily derive the local disper- 
sion relation for a single MHD disc and the local dispersion 
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Figure 24. The marginal stability curves for the case m = 0, 
/3 = 1/4, 5 = 1, a = 5. For both F = 0.6 and F = 1, the 
two cases of g = and g = 1.5 are computed. Compared with a 
smaller /3, the ring fragmentation regions rise up while the MHD 
Jeans collapse regions shrink. 



relation for a composite system of gravitationally coupled 
discs in the WKBJ regime can be expressed in a parallel 
form as equation 1981 . A s disc ussed for a single scale-free 
MHD disc bv lShen et all J2005I) . the WKBJ approximation 
is well justified for short wavelength perturbations (i.e., large 
£) but not as good for long wavelength perturbations. 

We comment on the MHD ring fragmentation instabil- 
ity here. In fact, this is the MHD generalization of the well- 
known Too mre instability characterized b y the so-called Q 
parameter ( Sa fronovll 19601 : lToomrelll964r) . Besides in con- 
texts of single discs, there have been numerous studies to 
identify such an effective Q parameter for a composite disc 
system felmegreenlll995l |jo3ll996t iLou fc Fanlll998bll an d 
in a singl e copl anar magnetized SID JLou fc Fanlll998al) . 
IShu" et alJ 1120001) noted that in a single SID, the minimum 
of the ring fragmentation branch is effectively related t o the 
Toom re Q parameter. For a c oplanar magnetized SID. iLoul 
l|2002ft and lLou fc Fan! i2002T) emphasized that the MHD 
ring fragmentation branch is also closely related to the gen- 



eralized Qm parameter jLou fc Fanll998"aF) . The conceptual 
connection of two instabilities com es from their depen dence 
on radial wavenumber. Physically, IShen fc Loul (2003) sug- 
gested a straightforward D criterion to determine the ax- 
isymmetric stability against ring fragmentation. In our com- 
posite scale- free discs with a coplanar magnetic field, we 
therefore interpret this MHD ring fragmentation as an ex- 
tension of the Toomre instability, with D 3 parameter being 
used for the stability criterion. In addition, this D s param- 
eter describes a global MHD instability criterion, while the 
Q parameter is only defined locally. 

When C2 1 = occurs for £ = £ c , one branch of the sta- 
bility curves diverges at £ = £ c . For a specified (3 value, the 
divergent point £ c depends only on F parameter contained 
in the coefficient C. 

For hydrodynamic pert urbations in c oupled scale-free 
discs without magnetic field. IShen fc Loul (|2004b) have dis- 
cussed the results thoroughly and their figures clearly show 
the influence of the two parameters a and 8 (see their Fig. 4 
to Fig. 8), for example, the increase of a or 5 lowers both the 
collapse and ring fragmentation regions. In other words, the 
increase of a or 8 reduces the danger of the Jeans collapse 
instability but makes the composite system more vulnera- 
ble to the ring fragmentation instability. By our exploration 
in the presence of magnetic field, this trend of variation re- 
mains. For this reason, we shall not repeat the analysis in 
this regard similar to theirs. We mainly focus on the case 
of 8 = 1, a — 5 and show the marginal stability curves 
fo r q = and q = 1.5 in one figure. Similar to the study 
of IShen et alJ i2005T) where global coplanar MHD pertur- 
bation structures in a single scale-free gaseous disc with a 
coplanar magnetic field are analyzed, we arrive at the same 
conclusion in our model of a composite system of coupled 
scale-free discs. The enhancement of a coplanar magnetic 
field reduces the danger of ring fragmentation in an appar- 
ent manner; and at the same time, it only suppresses Jeans 
collapse instabilities for the (3 range of— 1/4 < f3 < 1/4 
and tends to aggravate Jeans collapse instabilities for the /3 
range of 1/4 < /3 < 1/2. As has been already mentioned in 
IShen et alJ i2005t) . this collapse feature can be understood 
from the coupling between the surface mass density and the 
coplanar magnetic field of the background rotational MHD 
equilibrium. 

We also examine global MHD perturbation structures in 
a composite system of partial scale-free discs with the poten- 
tial factor F = 0.6. It is clear from Fig.gHland Fig. E^Elthat 
the presence of an axisymmetric dark matter halo suppresses 
both the ring fragmentation and Jeans collapse instabilities. 
This results from the fact that the dark matter halo exists as 
a stationary massive gravitational background. It can reduce 
the ring fragmentation because its gravitational force acts 
as centripetal force to resist such danger. It can suppress 
Jeans collapse because its stationar y background supports 
the stationary structure of the discs iMi ller et alJl97(j:l Hoh]| 
ll97ltlOstriker fc Peebleslll973l : Imiinev fc Tremainelll987tl 



5 SUMMARY AND DISCUSSION 

From several perspectives, this paper rep resents a general- 
ization and extens io n of pr evious w or ks ISver fc Tremaine 
(1996); IShu et alJ j2000h : iKalnaisl dl97ll) : lLoun2002); 
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i LoufcFanl (120021) : lLou fc Wul <l2005h : lLou & Zoul <2004 
l2006l) : IShen fc Loul <2004bl) : IShen et all <2005ty . We consider 
a composite disc system consisting of gravitationally cou- 
pled stellar and gas scale-free discs with a coplanar mag- 
netic field and construct global stationary coplanar MHD 
perturbation structures within such a composite disc sys- 
tem. More specifically, such discs are treated as razor-thin 
and scale-free so that the entire problem is two-dimensional 
with all background quantities varying in power laws of ra- 
dius r. For the background rotational MHD equilibrium, the 
surface mass density scales as Eo oc r~ a , the rotation curve 
scales as no oc and the azimuthal magnetic field scales as 
Bo on r * . The barotropic equation of state is expressed as 
II = ATE™. By the scale-free condition, the indices a, f3, 7, n 
are all expressed in terms of f3 by equation 1191 . The two 
discs are gravitationally coupled through the Poisson inte- 
gral and the scale-free condition requires that the physi- 
cal quantities of the two discs should share the same ra- 
dial scalings with different proportional coefficients allowed. 
Meanwhile, the potential factor F for a composite partial 
disc system is included to represent the ratio of the gravi- 
tational potential arising from the composite disc system to 
that arising from the entire system including a massive ax- 
isymmetric dark matter halo mass distribution which is as- 
sumed unresponsive to perturbations in the composite disc 
system. For convenience and clarity, several dimensionless 
constant parameters are introduced. The first one is the ef- 
fective rotational Mach number D (we omit either subscript 
or superscript 'g' or 's' for simplicity here) defined as the 
ratio of the disc rotation speed to the disc sound speed yet 
with an additional scaling factor (1 + 2/3)" 1 / 2 . In essence, 
this is a measure for the disc rotational velocity. The sec- 
ond one is q parameter corresponding to the ratio of the 
Alfven wave speed to the gas sound speed which provides a 
measure for the strength of the coplanar magnetic field. The 
other two parameters are 8 and r\ which represent the ratio 
of surface mass densities of the two discs and the square of 
the sound speed ratio (the velocity dispersion in the stellar 
disc mimics the sound speed in our formulation) of the two 
discs, respectively. For the convenience of figure presenta- 
tion, we also define a = I/77. These parameters completely 
characterize the behaviour of a composite disc system for 
the background equilibrium and are independent of each 
other; to a greater extent, they can be freely chosen but 
within a certain bound. In our model, there is a wide range 
of freedom to choose these parameters. Our goal is to study 
stationary coplanar MHD perturbations to construct global 
aligned and unaligned configurations described by the rele- 
vant dispersion relations containing all these parameters and 
to analyze partially the linear instability problem of such a 
composite disc system. 

Our model is highly idealized in several aspects, nev- 
ertheless, these simplifications are sensible and useful in 
making the problem tractable. First, we adopt the razor- 
thin approximation for the composite disc system so that 
vertical vari a tions across the d is c syst e m are ignor e d [e.g ., 
Qianl <1992l): ISyer fc Tremaind <1996ft : iFan fc Loul <1997l) : 
Shu et alJ (l2000D : lLou fc Fan! Jl998bF |. This is justifiable so 



dark matter halo. In reality, there exists another compo- 
nent of cosmic rays, important from the diagnostic perspec- 
tive for observing large-scale radio structures as revealed 
by synchrotron emissions from relativist ic electrons gyrat - 
ing around the galactic magnetic field (lLou fc Fan! 12003'). 
Technically, we need to deal with subsystems in a piece- wise 
manner and different subsystems will manifest different as- 
pects of the physical problem. In our model formalism, we 
mainly focus on the magnetized interstellar gases and the 
stars in a galactic disc system in the presence of an axisym- 
metric dark matter halo. This is typical for disc galaxies in- 
cluding our own Milky- Way. Thirdly, we adopt a simple fluid 
description for both the stellar and magnetized gas discs. In 
principle, a distribution function approach to the stellar disc 
would be more accurate and comprehensive. However, for 
a large-scale dynami cs without resonances etc., a fluid ap- 
proximation suffices (|Lin fc Shu||l964: Binnev & Tremainc 
[ 19871: iBertin fc Romedll98gl lLou fc SherJl2003l lLou fc Zoul 
2004; She n fc Loul l2004bl) . Fourthly, the basic problem is 
formulated here within the scale-free framework. Now all 
quantities are required to vary with power-law scalings in r, 
and the barotropi c equa t i on of state must be i m posed [e.g., 
ISver fc Tremaind (1199611 : IShen fc Loul <l2004bl) : IShen et alJ 



(2005)]. The scale-free condition further requires that the 
two coupled discs share the sa me radial scalings as w ell as 
the same barotropic index [e.g.. IShen fc Loul J2004b)]. An- 
other major constraint is that such scale-free configurations 



e.i 



; Lou fc Fan 



long as the radial scale of a galactic disc is much larger than 
the disc thickness. Secondly, we simplify a galactic system to 
just three major constituent components: a magnetized gas 
disc, a stellar disc (treated as a fluid) and an axisymmetric 



must b e stationary with oj 

|2002 |k lLou fc Shed diool: lLou fc Wul i2005ft: lLou fc Zou 
j2004l2006h : IShen fc Loul i2004blklShen et alJ (I2005M . Nev- 
ertheless, the model of coupled scale-free discs represents 
a significant extension of isothermal discs with /3 = 0. In 
particular, global MHD perturbation structures can be con- 
structed with the relevant MHD dispersion relations deter- 
mined analytically. Fifthly, the background magnetic field 
is presumed to be azimuthal and bears a radial power-law 
scaling r~ 7 . As the magnetic field distribution in disc galax- 
ies is generally complex, either the isopedic or the coplanar 
magnetic field geometry represents an idealization. Model 
analysis in terms of these two 'orthogonal' geometries pro- 
vide a theoretical basis for physical intuitions about mag- 
netic field effects. Our current model focuses on the case of 
a coplanar magnetic field, while for the case of a n isopedic 
magne t ic field, the intere sted reader is referred t o lShu et alJ 
<200CD : lLou fc Wul (120051) and lWu fc Loul(l2006h . 

Our model makes advances in the research in the MHD 
disc problem. By including two coupled scale-free discs with 
a coplanar magnetic field, our mo del contains variou s mode l 
ingredients of p r evious mod els of ISjgjLfcTremajrie] jl 996 ) ; 
IShu et alJ(l200d):lL^ll2002]):lLou fc Far] d2002[) : |Lou fc Zoul 
j2004l2006l) : IShen fc Loul <2004bl) : IShen et alJ ll200Eft . In a 
fluid-magnetofluid composite disc system with a coplanar 
magnetic field, we construct both aligned and spiral global 
MHD perturbation configurations and they are in accor- 
dance with previous works and meanwhile there are several 
new features in our model. With such an analysis, the MHD 
dispersion relations are derived analytically without any ap- 
proximation, meanwhile, the axisymmetric marginal stabil- 
ity criteria is also obtained analytically. Our model provides 
a new perspective for understanding large-scale MHD struc- 
tures in disc galaxies. We now summarize the main results 
in the following. 
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(i) Aligned Stationary MHD Configurations 
For the aligned cases, we derive the MHD dispersion 
relation and provide ample numerical examples concern- 
ing global MHD perturbation configurations and the corre- 
sponding phase relations among perturbation variables. In 
our nuid-magnetofluid model of two coupled scale-free discs, 
the MHD dispersion relation is a cubic equation in terms of 
y = Dg and there are three y = D 2 solution branches in 
general. By our analysis, the existence of the three branches 
is due to the gravitational coupling of perturbations in the 
two discs and the interaction of the magnetic field with the 
disc system [see the dispersion relation of equation 1651 1 
and they correspond to purely azimuthal propagation of su- 
per fast MHD density waves (sFMDWs), fast MHD density 
waves (FMDWs) and slow MHD density waves (SMDWs), 
respectively. Generally speaking, the yi branch (sFMDW) 
is always positive and the 1/3 branch (SMDW) is more often 
negative thus unphysical. Our interest lies mainly on the 
m = 2 (i.e., bar-like) case and we discuss the influence of 
several parameters. 

The magnetic field serves as a centripetal force when 
(3 < 1/4 (i.e., magnetic tension force is stronger than the 
magnetic pressure force) and as a centrifugal force when 
f3 > 1/4 (i.e., magnetic tension force is weaker than the 
magnetic pressure force) for the background MHD rotational 
equilibrium state. In the presence of coplanar MHD pertur- 
bations, there is a gravitational interaction between the two 
discs and there are also MHD interactions within the mag- 
netized gas disc. For cases of q 2 -C 1, the magnetic field 
has very weak influence on the composite disc system. The 
yz branch becomes t rivial. Our y \ and y 2 branches corre- 
spond to solutions of lShen fc Loul (l2004bl) . While for q > 1 
(i.e., the Alfven wave speed is comparable to or larger than 
the sound speed in the gas disc), the behaviour of solu- 
tion bra nches can be totally d ifferent from the unmagnetized 
cases of IShen fc Loul fe004bl) and magnetic field effect can 
become dominant, especially for small o (e.g., o ~ 1). In the 
presence of magnetic field, we observe a new feature that as 
(3 and q become large enough, there appears a 'convergent 
point' where for different 8 values some solution branches 
converge and then their relative positions are changed. Un- 
der the assumption that 77 ^ 1 (or a ^ 1, the gas sound 
speed is less than the velocity dispersion of the stellar disc) 
and when /3 < 1/4, we have D s < D g according to equa- 
tion (1221 : the magnetized gas disc rotates faster and the 
existence of magnetic field enlarges the rotational speed dif- 
ference between the gas and stellar discs. Yet for a relatively 
strong magnetic field and for f3 > 1/4, there are solutions 
with D a > D g and the stellar disc rotates faster instead (see 
Fig. Hi. 

In the aligned cases, the phase relations among pertur- 
bation variables are relatively simple because all quantities 
are real and they must be either in phase or out of phase. 
For the t/i branch, the perturbation mass densities in the two 
discs are always out of phase while for the other two branches 
2/2 and 2/3, the perturbation mass densities are generally in 
phase. However for the 2/2 branch with a strong magnetic 
field, the perturbation mass densities can also become out 
of phase. For the radial component of the magnetic field per- 
turbation, ii? and S 9 are always in phase for the 2/1 branch; 
they are in phase for the 2/2 branch and out of phase for 
the 2/3 branch in many cases but on the whole this phase 



relation depend on both q and o parameters (see Fig. HOB . 
Since Z = (4/3 — l)ii?/2m, the phase relation between gas 
mass density perturbation and the azimuthal component of 
the magnetic field perturbation is similar to that of S 9 and 
iR when j3 > 1/4 and has an additional phase shift n for 
(3 < 1/4. When [3 = 1/4 there is no azimuthal magnetic 
field perturbation. 

We also note that when f3 becomes sufficiently small 
(e.g., j3 = —0.24) with a mild magnetic field strength, there 
exists only one solution branch and the other two are com- 
plex conjugates. 

Regarding m dependence, there are not many new fea- 
tures except for m = 1. In this case, the solutions are quite 
different from m = 2 modes. In particular, it is now possible 
that the solution branches intersect and also there can be 
no physical D 2 solutions. 

(ii) Stationary Configurations of Logarithmic Spirals 

The MHD dispersion relation for the spiral cases is 
quite similar to that for the aligned cases in form. Differ- 
ences only occur by a systematic replacement of notations. 
Therefore, there are many common features in solution be- 
haviours. Physically, logarithmic spirals involve both radial 
and azimuthal propagations of MHD density waves, while 
the aligned cases inv olve purel y azimuthal propagations of 
MHD density waves (lLoull2002f) . 

In the m — 2 case, the 1/3 branch is due to the exis- 
tence of magnetic field and seems to bear little relation with 
the spiral pattern (i.e., it almost does not vary much for 
arbitrary £). This SMDW is in most cases negative thus un- 
physical. Those physical SMDWs in the spiral cases appear 
when (3 is not too small (i.e., j3 ^ 0) and a is small enough 
(i.e., g = 1). The 2/1 and 2/2 solution branches are monoton- 
ically increasing with radial wavenumber £, suggesting that 
to maintain a stationary tight-winding logarithmic spiral 
pattern, a fast disc rotation speed is required for sFMDWs 
or FMDWs. The increase of the magnetic field strength gen- 
erally raises the solution branches, while a large o tends to 
lower the solution branches. The 'convergent point' appears 
more often in the logarithmic spiral cases as £ varies. 

We now describe the phase relations among perturba- 
tion variables for logarithmic spirals. The phase relation be- 
tween perturbation surface mass densities remain either in 
phase or out of phase, and the behaviour is similar to those 
in the aligned cases, but the phase relation between surface 
mass density perturbation and the magnetic field pertur- 
bation becomes fairly involved. The dependence of S 9 /(iR) 
and S 9 /Z on £ is described by a complex expression. Now 
the azimuthal magnetic field perturbation is Z — £LR/m 
which is quite different from that of the aligned cases and 
is dependent on the sign of £ (i.e., for leading waves and 
trailing waves the azimuthal magnetic perturbation have a 
phase difference of 7r). In the case of £ = for 2/1 and most 
of 2/2 branches, the ratio S 3 /(iR) is positive, i.e., the ra- 
dial magnetic field perturbation lags behind the gas surface 
mass density perturbation by n/2 in phase (no azimuthal 
magnetic perturbation). As £ increases, the imaginary part 
of this ratio decreases in general while its real part does not 
vary significantly, which indicates that the phase of the ra- 
dial magnetic field perturbation is catching up with the gas 
surface mass density perturbation (meanwhile the azimuthal 
magnetic perturbation becomes ahead of S 9 from the orig- 
inal in-phase relation). In the tight- winding approximation 
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(i.e., £ — > oo), R and S 9 are in phase and the phase of az- 
imuthal magnetic perturbation tends to be ~ tt/2 ahead of 
S 9 . We also find that in the /3 = 1/4 case, for a small a = 1 
and a large <? = 1.5, there is an abrupt phase shift of n in j/2 
branch as £ increases from 0. In the physical j/3 branch for 
SMDWs, the phase relation between S 9 and R or S 9 and Z 
has an approximate difference of it compared with that of 
the yi branch. 

(iii) Axisymmetric Marginal Stabilities and Magnetoro- 
tational Instabilities 

We explore the unaligned m = perturbations referred 
to as the marginal stability curves where two regimes corre- 
sponding to the Jeans collapse instability and the ring frag- 
mentation instability appear. The criterion for the global 
ring fragmentation instability is an MHD generalization of 
the local Toomre Q-criterion. 

In our fluid-magnetofluid model, we arrive at the same 
conclusions as discussed by Shen & Lou (2004b) and Shen 
et al. (2005) that large a and large 8 lower both the collapse 
and ring fragmentation regions. In other words, the increase 
of o~ or S reduces the danger of Jeans collapse instability 
but makes the system more vulnerable to ring fragmenta- 
tion instability. The enhancement of magnetic field reduces 
the danger of ring fragmentation, suppresses Jeans collapse 
instabilities for —1/4 < /? < 1/4, and aggravates Jeans col- 
lapse instabilities for 1/4 < /3 < 1/2. 

An interesting result is that when we examine the m — 
case for both aligned and unaligned perturbations, we find 
that the aligned m — case is just the limit of the unaligned 
m — case with the radial wavenumber £ — > (i.e., the 
breathing mode) which does not merely represent a rescaling 
of the background MHD rotational equilibrium. 

We discuss briefly on the well-known magneto- 
rotational instabilities (MRIs) in a weakly magnetized disc 
system. In general, for a disc of finite thickness, if its angu- 
lar velocity decreases outward with a weak magnetic field, 
the disc rotational profile becomes unstable caused by MRIs. 
Such MRIs are widely invoked to generate and sustain MHD 
turbulence in the disc and result in an outward transport 
of angular momentum. Both local and global MRIs appear 
in extensive numerical simulations. Although MRIs are ini- 
tially explored with a vertical magnetic field, it can also 
manifest when the background magnetic field is azimuthal 
l)Balbus fc HawleJl998UOgifvie fc Pringlelll996h . Although 
a purely azimuthal magnetic field is less effective in produc- 
ing an MHD turbulence, numerical simulations do show that 
a weak vertical magnetic field can be produced from an az- 
imuthal magnetic field and the entire system then behaves 
simila r to a system with a purely vertica l magnetic field 
ijHawlev et al.ll99,4lBalbus fc Hawlevll998h . Such MRIs are 
suppressed when the magnetic field becomes so strong such 
that magnetic a nd rotational energies ap proach an equipar- 
tition condition IjQgilvie fc PringleMl996r) . Since only razor- 
thin discs are considered in our current formulation, MRIs 
do not appear in our model. Nevertheless, in reality, the fi- 
nite thickness of a gas disc can induce MRIs and result in an 
outward transport of angular momentum. Moreover, locally 
produced MRIs in a galactic disc may further complicate 
the appearance of large-scale spiral patterns discussed here. 

(iv) Partial Discs with a Dark Matter Halo 

We apply previous analysis to a composite system of 
partial discs with F < 1. For the m = 1 case, the po- 



tential ratio parameter F plays an important role to affect 
the behaviour of Dg solutions. In the aligned m = 1 cases, 
when /3 is not too small (e.g., (5 > 0), there is a 'divergent 
point' F c determined by /3 where the y\ branch diverges and 
physical solutions appear only for F ^ F c . When [3 is very 
small (e.g., (3 — —0.24) together with other parameters (i.e., 
a — 1.2, q — 0.4), the only physical solutions lie in the 
interval < F ^ F C 2 [see panel (b) of Fig. 1111 . For station- 
ary logarithmic spirals with m = 1, the value F c depends 
on both P and £ and we would propose that spiral galaxies 
with stationary patterns of one arm and a large proportion 
of dark matter halo should have large radial wavenumber 
(i.e., tightly wound one-arm spiral). 

When m ^ 2 in both aligned and unaligned cases, the 
dependence of D$ solutions on the potential parameter F is 
generally weak. Nevertheless, F may change phase relations 
of perturbation variables significantly. First, we discuss the 
aligned cases. When /3 is small (e.g., j3 = —0.24), the ratio 
S 9 / S" for the 1/1 branch can be positive for small £ and small 
a although for a full disc this ratio remains always negative. 
For the phase relation between S 9 and R, one feature is that 
when q is large, R can be zero [i.e., a singularity in the ratio 
S 9 /(iR)] for a certain £ c . Meanwhile, there are also some 
variations in the phase difference for partial discs which are 
clearly shown in panel (c) of Fig. 1191 For j3 = 1/4 case and 
in full discs, the 1/2 branch has negative S 9 / S s ratio when 
a = 1. In a partial composite disc system with F — 0.5, for 
the entire 1/2 branch the ratio S 9 / S" is positive. Meanwhile, 
for both yi and j/3 branches, the magnitude ratio \S 9 /S S \ 
is much larger than that of full disc cases. For the phase 
relation between S 9 and R, we no longer have R = 0, but the 
phase variation is still evident (see Fig. l20l for a comparison). 

Finally, in the m = case with purely radial oscilla- 
tions in the composite MHD system, the presence of a dark 
matter halo mass greatly reduces the danger of both Jeans 
collapse and ring fragmentation instabilities. In other words, 
a composite system of partial discs appear more stable than 
full discs. 
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APPENDIX A: A FURTHER DISCUSSION FOR 
THE ALIGNED AXISYMMETRIC CASE 

As already noted in the main text, the m = case should 
be calculated in the process of first setting m = and then 
taking the limit u — > 0. We therefore start from equations 
lUTfl l - (O for our calculation. Now with S, E oc r~ 1_2/3 
in equations I49H and H53L we have U tx r; it then follows 
from equations (15 IP and 1531 1 that J oc r 1- ' 3 and *P oc r~ 2/3 
as before. Under such circumstances, the radial scalings of 
certain perturbation variables are different from before and 
the scale-free condition is invalid if uj 7^ holds. However, 
we find that in the limit u> — > 0, we must require U — > 0. 
Now when the first term on the left-hand side of equation 
<50> together with the first term on the left-hand side of the 
second relation of equations I53H approach zero, the scale- 
free condition becomes satisfied again. Meanwhile, U/uo 7^ 
in such a limiting process, leading to a solution for this case. 
We show the result of this limiting process below. 



^-fltfrV = 2/3E*** 



(1 - 2/3) 

2(1 -mlr 2 
(1-2/3) ' 



Gj(4/3-l)/3 
(1 - 2/3) 



5' ! ' 



(Al) 



2/3Eo* s 



By a substitution of equations I36H and 1541 into equation 
JAll l and a further simplification, we obtain 



2 /3 a 2 + + ° aI3{1 ~ 4P) 



(1 - 2/3) 



27rGr(2/3P )Sg 



S 9 - 2-KGr(2l3V )T, 9 S s = , 



(A2) 



Wa 2 s + 2(1 (1 - 2irGr(2pP )E° 

- 2TxGr{2f3Vo)T, s S 3 = . 



The two relations of equation 1A21 cannot be simultaneously 
satisfied unless the coefficient determinant of equation \K2\ 
vanishes; for non-trivial solutions, this requirement leads to 
the dispersion relation in the form of 



[V 9 + T - 27rGr(2/3P )Sg][V s - 27rGr(2/3;Po)£o] 
= 47r 2 G 2 r 2 (2/3Po) 2 SgEg , 

where, for notational simplicity, we define 

V = 213a 2 + 2(1 - /3)fi 2 r 2 /(l - 2/3) , 
^=Cl/3(l-4/3 2 )/(l-2/3) . 



(A3) 



(A4) 
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Substitutions of equations 1201 — 123H in equation 1A3II yield 



2/3 F(D 2 g + 1 + Q)6 2(1 - (3)D 2 g 



2/3Q 



1+2/3 



1 + 6 



1 - 2/3 



2/3 



2/3 



F(D 2 + 1) + 2(1 - /3)£ 2 



1 + 5 



1-2/3 



(A5) 



(£> 2 + l)(£> 2 + l + Q)F 2 a 
(l + <5) 2 



where the magnetic field effect is contained in the following 
Q parameter 



Q = (4/3 - l)g 2 /(4/3 + 2) 



(A6) 



By evaluating equation IA5t using radial force balance rela- 
tion [1221 . we obtain a quadratic equation of D 2 which has 
two roots. To show the physical meaning of this equation, we 
intend to compare this equation with equation I99H where 
axisymmetric marginal stabilities are discussed. We expect 
that by taking £ = and j3 — 1/4, equations lA5t and 
19911 should coincide and this is indeed the case by noting 
AT m (l/4, 0) = Pm(l/4), Q(l/4) = 0, Ai'(0) = 1/4. This 
result unambiguously shows that the aligned m = case 
does actually correspond to the limiting case £ — > of the 
m — logarithmic spiral case. Since the potential-density 
pairs carry different scales, only in the case /3 = 1/4 do they 
coincide. 1 Now it turns out that the aligned m = case is an 
extension of the marginal stability curves in the limit £ — > 0. 

This result may be understood as follows. In the aligned 
m = mode, the equation of continuity is invalid unless 
in the limit of uj = 0. Since £i(r) = er~ 1 ~ 2/3 exp(iujt), the 
surface mass density of the whole disc varies synchronously 
if tv 7^ 0. However, we may regard such mode as the limiting 
case of the spiral m = case with radial wave number £ — > 
and write 



Ei(r) = er~ 1_2,3 exp(£lnr + iwt) 



(A7) 



Now when £ and u is nonzero but extremely small, the equa- 
tion of continuity still holds. Here, the surface mass density 
perturbation scale differs from that of the marginal stability 
curves in Section 4.3 unless /3 = 1/4. Thus in this perspec- 
tive, the result may be regarded as the marginal stability 
relation with perturbation scale Ei oc r ~ 1 ~ 2 i 3 + 1 ^ for small £ 
(i.e., widely open spirals) cases and is an extension of our 
discussion on marginal stability problem in Section 4.3. 



APPENDIX B: CALCULATIONS 

OF DISPERSION RELATION 
FOR THE ALIGNED CASES 

We evaluate the MHD dispersion relation by computing 
the coefficient determinant of equations 15611 . 1571 . 1581 
and 1601 1 for a specified background equilibrium described 
by equations 121H — 1231 1. The explicit coplanar MHD per- 
turbation equations in terms of characteristic parameters 



1 In potential-density pair 1761 for logarithmic spirals, the per- 
turbed surface mass density does not scale the same as the back- 
ground surface mass density unless /3 = 1/4. 



0, q, m, 5, t], D g , 
step, we obtain 



F are summarized below. In the first 



mD g S 3 
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2m(l + 2/3)D g r 1 +f } 



M^J 9 -4^G/3P m (/3)S 



(B2) 



[2/3P (/3)][4/3 + (l-4/3)g 2 ](l + a) 
(1 + 2/3) \D 2 + 1 + (4/3 - l)g 2 /(4/3 + 2)] 



S 9 =0 



(1-/3) 
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r i+/3 (4/3 + 2)D 9 r 1 +f 
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[D 2 + l + (4/3 - l)g 2 /(4/3 + 2)]F6 

for the magnetized gas disc and 

[m 2 - 2(1 - (3)]D 2 S 1 
[m 2 + 4/3(1-/3)] (1 + 2/3) 



it/ 9 

2-KGm^P m (P) 
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+ 



(D 2 + l)Vm(P) 



S 9 + 



(D 2 B +l)V m (P) 



(B4) 



S s 







[2pP (p)](l + 5)J [2/3P (/3)](l+«5) 

for the stellar disc. In addition to the notations already de- 
fined by equations 16611 and 1A6II . we here introduce 

C(/3) = 2/3/(1 + 2/3) , 



H(D t ,p,q) = (D 2 + 1 + Q)F/P (J3) , 
and note from equation 122|l that 



(B5) 



D 



vHPo - 1 • 

Using these newly defined notations, equations 1B11 — 1B4> 
take the forms of 



_ SAgHSiU 9 
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r 2+/3 

for the magnetized gas disc and 



+ 



C(l + S) 
yHP m ((3) , 
G(l + <5) ' 
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34 Y.-Q. Lou, X. N. Bai 



for the stellar disc. We rearrange the terms in these four 
equations in the order of (S 3 , ill 3 , J 3 , S s ) to identify and 
manipulate the coefficient determinant by drawing several 
factors and by properly adding and subtracting lines and 
rows. The simplified determinant appears as 



(BIO) 



an 


ai2 


dl3 





bn 


bi2 
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bi4 


Cll 


(12 
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Cl4 


du 








du 



where the eleven coefficients are given explicitly as 



an 
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3SD 2 n 
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6n = -H/3Vm + 13(1 + 5)(C - Q)/S , 

b 12 = m 2 D 2 a + W- 1 / : 2m " Q , 6i4 = -H(3Vm , 



an 
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ci4 = -HP m , du = v HV m A m /[C(l + 6)] , 
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where the eleven coefficients are explicitly defined by 
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C (1 + o) 

It is then straightforward to evaluate this determinant nu- 
merically and the results are shown in Section 4 for loga- 
rithmic spirals of MHD perturbations. 

This paper has been typeset from a TjjjX/ I^TjtX file prepared 
by the author. 



With these definitions, it is straightforward to evaluate the 
determinant numerically and the result is shown in Section 
3. We note that in acquiring the coefficients of the polyno- 
mial equation of y in the form of equations 1681 and 16911 . it 
is somewhat onerous and tedious to combine fragment terms 
(e.g., m 2 , /3 2 , f3) to form extra coefficients Am and B m etc. 



APPENDIX C: CALCULATIONS 

OF DISPERSION RELATION 
FOR LOGARITHMIC SPIRALS 

Since the MHD perturbation equations for logarithmic spiral 
perturbations parallel those for aligned perturbations, the 
procedure of derivations is strikingly similar. We will not 
repeat the steps but directly summarize the key results of 
our calculations. 

Similar to the aligned case, the main purpose is to ex- 
press the coefficient determinant of equations 15611 . 15711 . 
15811 and 16011 in reference to the background equilibrium 
described by equations 1211 1 — 12311 in terms of the dimen- 
sionless parameters ((3, q, m, S, r), D g , F). Following the 
same steps of Appendix B and using the same notations, we 
show the simplified determinant below. 
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